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Abstract 



Analytic Approaches to the Study of 
Small Scale Structure on Cosmic String Networks 

by 

Jorge V. Rocha 



We present an analytic model specifically designed to address the long standing 
issue of small scale structure on cosmic string networks. The model is derived 
from the microscopic string equations, together with a few motivated assumptions. 
The resulting form of the correlation between two points on a string is exploited 
to study smoothing by gravitational radiation, loop formation and lensing by 
cosmic strings. In addition, the properties of the small loop population and the 
possibility of detecting gravitational waves generated by their lowest harmonics 
are investigated. Whenever possible, we compare the predictions of the model to 
the most recent numerical simulations of cosmic string networks. 
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Chapter 1 
Introduction 



The subject of cosmic strings lies in the interface of many branches of Physics: 
first proposed in the context of field theories and in close analogy with condensed 
matter systems, their existence would be relevant for cosmology and imply de- 
tectable gravitational effects. More recently the connection with string theory has 
also been made, which resulted in the renewal of interest on the field. 

1.1 Cosmic strings: what are they? 

The idea that topological defects could form during a cosmological phase tran- 
sition was first suggested by Kibble in 1976 [1]. The picture behind this proposal 
was that of a gauge theory that undergoes a symmetry breaking transition at some 
critical temperature. In a hot Big Bang model, the temperature of the universe 
grows without bound as we go back in time. At very high temperatures the gauge 
theory finds itself in a symmetry preserving configuration, but as the universe 
expands and consequently cools down, the critical temperature is reached and the 
theory enters a broken symmetry phase. 

In particular, whenever the field theory has a spontaneously broken U{1) sym- 
metry there are classical solutions which are localized in two spatial dimensions. 
Such solitons are present for example in the Abelian Higgs model, where they are 
known as Nielsen-Olesen vortices [2]. This model couples an Abelian gauge field 
to a complex scalar field with a typical "Mexican hat" potential for which there 
is a ring of degenerate vacua (see Fig. II. ip . The scalar field configuration for the 
vortex maps circles in space (with sufficiently large radius) to the ring of the po- 
tential minima. Continuity implies that somewhere in the middle the field sits at 
the local maximum of the potential. Analogously, if we lay down a blanket in such 
a way that its rim encircles, say, a pillow, then some portion of the blanket must 
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Figure 1.1: The usual quartic "Mexican hat" potential for the complex Higgs field 
0(x) is shown on the left. A vortex string configuration maps spacetime points 
belonging to the core of the string (on the right side of the figure) to the local 
maximum of the potential V{(j)) and large circles surrounding the string to the 
ring of degenerate minima of V. 

cover the pillow. So the solution has a core where the energy density is localized. 

One can now easily imagine extending this solution in one more dimension, 
thus obtaining a string. The discussion above makes it evident that these defects 
are classified by their winding number, the (integer) number of times the field 
wraps around the ring of vacua as we trace out a closed path in physical space. 
The stability (at least in the absence of sources) of these objects stems from this 
topological nature. Summarizing, cosmic strings are linear threads of energy with 
cosmic-scale extensions that are produced in cosmological phase transitions. Such 
objects can be very massive as their core retains the typical energy density of 
the universe before the phase transition. Thus, cosmic strings represent remnants 
from the early universe. 

Of course, other topological defects are possible. The simplest one is the do- 
main wall obtained as the 'kink' solution for a real scalar field theory interpolating 
between the two minima of a potential with discrete Z2 symmetry, trivially ex- 
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tended in the two transverse directions. Monopoles, which are point-hke defects, 
also arise in theories for which the manifold of degenerate vacua is topologically 
a sphere. All the possibilities are nicely classified by homotopy groups. 

The crucial point is that these defects inevitably form in cosmological con- 
texts by the Kibble mechanism. As the temperature drops below the critical 
value, regions of space-time in the broken symmetry phase start nucleating. The 
choice of vacuum by each region must be uncorrelated on distances larger than 
the cosmological horizon. Typically, these nucleated bubbles expand with veloci- 
ties approaching the speed of light and eventually start smashing into each other. 
When that happens, occasionally strings are formed when the energy density in 
its core gets trapped by the non-trivial winding of the field around it. The con- 
servation of this winding number implies that these strings cannot end as this 
would require an excursion of the field over the top of the potential and such large 
fluctuations are prohibited once the temperature has dropped below the Ginzburg 
temperature. So they either have infinite extension or they form closed loops. 

The same rationale applies to defects of different dimensionality but domain 
walls and monopoles are both cosmologically dangerous: their energy density 
would dominate over the driving source of space-time expansion during part of the 
history of the universe, therefore ruining the cosmological evolution. To prevent 
this catastrophe we postulate a period of accelerated expansion of the universe, 
known as inflation, thereby diluting all defects formed earlier. Hence, cosmic 
strings will only be relevant if they form after or at the end of inflation. 

1.2 Cosmic strings in our universe? 

Following the seminal paper by Kibble there was a period during which the 
field was dormant but four years later Zel'dovich |3j and Vilenkin [4] brought an 
additional thrust to the subject, by pointing out that cosmic strings produced 
during a phase transition at a Grand Unification Theory (GUT) energy scale rjsB 
(~ lO^^GeV) would give rise to density fluctuations of the right magnitude to ex- 
plain galaxy formation. Indeed, the tension fi of the string is proportional to the 
square of the symmetry breaking scale. Multiplying by Newton's constant G and 
the appropriate power of the speed of light c we obtain the all-important dimen- 
sionless quantity Gfi/c^ ~ rfg^/Mpi. Here Mpi ^ 1.3 x lO^^GeV/c^ is the Planck 
mass, representing the scale at which gravitational and quantum effects become 
equally important. One then concludes that GUT scale strings would generate 
effects of order G/z/c^ ~ 10~^. This combination is proportional to the fractional 
change in density that would be generated by cosmic strings and the value just 
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given roughly corresponds to what is needed to seed galaxy formation. Further- 
more, it has been shown that the production of cosmic strings in supersymmetric 
GUT models is quite generic [5J. 

Much work has been done since then in order to study properties and conse- 
quences of cosmic strings in the universe0 However, the interest in the subject 
started to fade away by the end of the nineties with the advent of experiments 
devoted to the study of the anisotropics in the Cosmic Microwave Background 
(CMB), namely COBE (Cosmic Background Explorer) and more recently WMAP 
(Wilkinson Microwave Anisotropy Probe). The decline was partly caused by the 
very good agreement between the angular power spectrum of the CMB revealed 
by the measurements and the predictions of the competing theory of inflation. 

In the inflationary paradigm, the density perturbations are generated by quan- 
tum fluctuations occurring during the period of accelerated expansion. Modes 
with wavelengths greater that the Hubble distance H{t), whose inverse gives the 
expansion rate of the universe, are frozen in. When traced back in time, comoving 
distance scales observable today (and therefore smaller than H) cross the Hubble 
radius and remain frozen for some time before they re-enter. Different scales stay 
frozen for different amounts of time in a very specific way and this in turn results 
in a distribution with characteristic peaks for the angular power spectrum of the 
density perturbations. 

On the contrary, in the cosmic string scenario the density anisotropy is a direct 
consequence of the inhomogeneity of a universe filled with a network of strings. 
The original proposal was that matter accreted around closed loops of cosmic 
strings. At that time it was thought that the loop number density was very close 
to the density of galaxies, suggesting a real connection between the two entities. 
However, simulations performed by the end of the eighties [H El [TOl |TT] showed 
that this was not the case: typically the string networks evolve into configurations 
with a much denser gas of loops, with sizes much shorter than the horizon scale 
and correspondingly increased velocities. So the picture was replaced by the idea 
that structure formed in the wakes of moving long cosmic strings. This is made 
possible by the gravitational properties of one-dimensional defects: a straight 
string produces a conical deficit in the geometry of the transverse directions so as 
it passes between two initially static objects it generates a relative attractive mo- 
tion. Still, this mechanism of generating density perturbations predicts a broad, 
nearly featureless, distribution for the CMB angular power spectrum and it be- 
came apparent that cosmic strings were disfavored already with the data collected 
by COBE [12]. 

^The interested reader is directed to Refs. [51 [7] for nice detailed expositions. 
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So cosmic strings were dismissed as the main source of density perturba- 
tions responsible for structure formation. Indeed, their contribution to the CMB 
anisotropy has been constrained to be no more then roughly 10% [13]. In a sense, 
it can be said that the inflationary theory had more chances of taking the lead as 
we have already mentioned that even in the case of a universe containing topologi- 
cal defects, we still need to appeal to inflation to avoid cosmological catastrophes. 
Nevertheless, cosmic strings are not ruled out and in fact recent studies [T^ [T5] 
indicate that a small contribution from these defects can actually provide better 
fits to the data, possibly explaining the excess of power in the anisotropy at small 
scales, for example. 

1.3 Cosmic superstrings 

More recently, the possibility of the emergence of cosmic strings from models 
of inflation in string, theory, namely the brane inflation scenario [16], has driven 
a revival of the fieldo In this incarnation, topological defects may arise as funda- 
mental strings (F-strings), Dl-branes (D-strings), bound states of F- and D-strings 
or higher dimensional branes wrapping various cycles in the extra dimensions [T^ . 

When the possibility of superstrings taking the role of cosmic strings was 
considered in 1985 by Witten [19] it did not seem promising: the fundamental 
strings, having tensions of order the Plank scale, would generate anisotropies in 
the CMB much larger than those observed. In addition, these objects appeared to 
be unstable: in open string theories a long string would simply break into many 
small strings and closed strings would collapse due to the force exerted by the 
tension of a domain wall bounded by the string itself. 

The realization that other (non-perturbative) objects exist in string theory 
and that compact extra dimensions can include throat regions in their geometry, 
which would suppress the tension of strings lying therein, provided a much more 
fertile ground for cosmic superstrings. Indeed, in brane-inflationary models it is 
the inter-brane distance that drives inflation as two branes approach each other. 
Eventually a tachyonic mode appears and inflation ends when the branes collide 
and annihilate. The initial system contains U(l) gauge flelds and this symmetry 
is broken at the end of the process. So again one expects formation of lower 
dimensional defects. The types of extended objects we are left with depend on 
the particular D-brane model but, interestingly, from our four dimensional point 
of view one can (and even must) obtain only string-like defects [20| [2n [22| [23] . 
Therefore, the creation of cosmic superstrings at the end of inflation is expected 

^Reference [T^ gives an excellent introduction to the subject of cosmic superstrings. 
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whereas domain walls and monopoles are simply absent. 

It should be noted that the properties and stability of these cosmic super- 
strings are model-dependent. Depending on the D-brane geometry that pro- 
duces inflation, the resulting string tension can be anywhere within the range 
10~^^ < Gyu/c^ < 10~^ [2H [22] . The issue of stability has been reexamined in 
the light of the modern understanding of string theory in Ref. [18] where it was 
found that many models can accommodate cosmic superstrings that are at least 
metastable. 

1.4 Networks of cosmic strings 

When a network of cosmic strings is formed it is expected to resemble a tangle 
of random walks on scales greater than the correlation length of the Higgs fleld. 
This is conflrmed by numerical simulations [2l]. Some strings form closed loops 
but in an inflnite universe strings of inflnite extension (often called long strings) 
will also exist. After formation the network will not remain static. First of all, 
the strings have tension and so the segments will start moving around. On the 
other hand, spacetime itself is evolving and so the network is initially in an out- 
of-equilibrium state. 

The existence of cosmic strings populating our universe would generate ef- 
fects we could detect, in principle. Besides influencing the CMB angular power 
spectrum as already mentioned, other effects are also predicted: the conical na- 
ture of spacetime around a string gives rise to very peculiar gravitational lens- 
ing properties [2S] and this should also affect the CMB temperature pattern 
by introducing step-like discontinuities [26]; their gravitational radiation prop- 
erties leads to strong bursts of gravitational waves originating from special events 
(cusps and kinks) along the strings [27]; also, the presence of a gravitational 
radiation stochastic background would distort the very precise time-periodicity 
of pulsars [2B]; flnally, particle emission from annihilating strings, in particular 
gamma ray bursts [29], is a possibility as well. 

The lack of any such observation until present date translates into constraints 
on the cosmic string properties, typically an upper bound on the dimensionless 
string tension G/z/c^. The most stringent current bounds come from matching the 
CMB power spectrum [20], which gives G^/c^ < 2.7 x 10^'' at 95% confldence level. 
The effect of cosmic strings on the angular power spectrum comes mostly from 
the long string population and depends primarily on the large scale properties 
of the network. The reliability of the above bound stems from the fact that 
the large distance structure of cosmic strings is well understood. A recent study 
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of the gravitational lensing effect reports the constraint Gfi/c^ < 3.0 x 10^'' at 
95% confidence level [31] and this limit should be easily improved by an order of 
magnitude in the near future [32]. Measurements of pulsar timing exclude values 
of the dimensionless string tension above 10~^ [33j but these methods depend on 
more debatable network properties, most notably the typical size of loops formed 
by self-intersections of long strings and their velocities. 

The previous paragraph illustrates the importance of understanding the evo- 
lution of cosmic string networks for establishing constraints on the parameters. 
Another aspect is that knowledge about these systems provides a better guide 
to search for cosmic strings in our universe. Finally, if one day such topological 
defects are discovered, we would hope to be able to interpret the data and obtain 
valuable information about the microscopic theory that supports those extended 
objects. For example, one would wish to know if the detected strings (or rather 
their effects) were field theoretic cosmic strings or string theory counterparts. 

Fortunately, there are some differentiating features. In particular, when two 
segments of field theoretic vortex strings intersect the outcome is deterministic 
and for most of the range of parameters, namely the relative velocity, it results 
in a reconnection. On the contrary, in string theory this intercommutation is a 
quantum process and the probability of reconnection can be highly suppressed. 
Also the presence of bound states in string theory permits the existence of Y- 
junctions in the string network and such features are not possible for ordinary 
cosmic strings. 

As it turns out, the evolution of cosmic string networks is a notoriously dif- 
ficult problem. The string equation of motion in flat spacetime is a simple wave 
equation but once we consider expanding spacetimes the equations of motion be- 
come non-linear. In addition, the evolution of these systems is also influenced 
by intercommutation events (which generates kinks along the strings [31]) and 
gravitational radiation (which shortens the strings). As we will discuss, gravita- 
tional radiation effects are naturally suppressed by (powers of) Gfi/c^. Given the 
current bounds on the dimensionless string tension one concludes that there are 
large ratios of length and time scales involved in the evolution of cosmic string 
networks, rendering a full numerical analysis impracticable. 

The large scale properties, i.e. characteristics of the network on scales compa- 
rable to the horizon distance, have been fairly well understood from the earliest 
simulations [TUl [HI E] . As long as there are a few long strings at formation of the 
network, the late time configuration approaches the so-called scaling regime, in 
which all length scales grow proportionally with time. Independently of the initial 
conditions there will be a few dozen long strings crossing any horizon volume and 
on these large scales the strings trace out Brownian trajectories. Given the im- 



7 



CHAPTER 1. INTRODUCTION 



portance of the scaling regime we dedicate the following subsection to describing 
its nature. 



1.4.1 The scaling regime 

The idea of a scaling regime for a cosmic string network was first introduced 
in [35]. In addition to its obvious attractiveness, it was confirmed by the early 
simulations [361 EH [lOl El El [H] • An important consequence of the scaling regime 
is that the energy density in strings remains a (small) fixed fraction of the back- 
ground energy density in both radiation and matter cosmological it red- 
shifts with the same power of the scale factor in each epoch. The by now consen- 
sual value for the contribution from long cosmic strings to the total energy density 
is roughly GOG^/c^ (lO^Gyu/c^) times the matter (radiation) energy density in a 
matter (radiation) dominated universe [TOl [8l \TT \ [37 1 l38j. 

Several scales can be defined for a cosmic string network. Among them, the 
characteristic length of the network plays an important part. It is defined as the 
length scale L such that a typical volume of the network contains a length L 
of long strings. Another length scale is the persistence length ^, i.e. the distance 
along the string beyond which points are uncorrelated. In any case, we say that 
a quantity (with units of length) is scaling if it remains constant in units of the 
cosmological time t as the network evolves. Consequently, in a scaling regime the 
evolution is self-similar, in the sense that the network at a given time will resemble 
itself at a previous time after all lengths are scaled down appropriately. 

Thus, in a scaling regime we have L = jct, for some constant 7, and the energy 
density in long strings becomes 

Poo = ^ = ^ • (1.1) 

The consensual value for the proportionality constant is 7"^ ~ 3 in a matter- 
dominated universe and 7"^ ~ 10 in a radiation-dominated universe [SI [HI [33 
[38] . Similarly, the persistence length ^ is also observed to approach a scaling 
regime [391 EZ] • 

A crucial feature of cosmic string networks is that the scaling regime is an 
attractor. The intuitive picture behind this is the following. If the long strings 
simply grew with the expansion of the universe their density would decrease as 
a(t)~^, where a{t) oc t" is the scale factor and z/ = 1/2 for a radiation-dominated 
universe and i/ = 2/3 in a matter-dominated era. Note that such a density is 
growing with respect to a scaling density, which would decay as If we start 
out with too many long strings the probability of an encounter between two of 
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them will be large, leading to a bigger loss of long strings into loops. On the other 
hand, if we begin with a sparse density of infinite strings intercommutations will 
be rarer events and the long string density will increase relative to the scaling value 
and eventually reach it. Therefore, the late time scaling behavior is independent 
of the initial conditions as long as there are a few long strings to begin with§ 

More recently the issue of scaling in cosmic superstring networks has also been 
studied numerically [HI HJl HQl US]. The existence of Y-junctions could possibly 
lead the network to freeze but the scaling regime seems to be a robust prediction 
in this case as well. In fact, there are indications that the late time evolution of 
cosmic string and superstring networks might be very similar [H]. If this is the 
case, the results presented in this dissertation should also be applicable in the 
cosmic superstring scenario. 

1.4.2 Small scale structure 

By now, the evolution of structure on cosmic string networks on large scales 
is a well understood problem. The situation gets more complicated and less clear 
once we turn to the short distances. Indeed, the early simulations [IDl El flT] 
showed evidence for the build-up of small scale structure along the strings, caused 
by kinks generated during intercommutations, and for the presence of a gas of 
tiny loops with sizes accumulating at the smallest possible scale - the resolution 
scale. To date, there have been many analytic and numerical studies of cosmic 
string networks, but no single work has yielded completely satisfactory results. 
As mentioned above, the large ratios of length and time scales involved make a 
full numerical treatment impracticable. On the other hand, analytic methods are 
difficult because of the highly non-linear nature of the system. The hope is that 
a careful combination of analytic and numerical approaches will lead to a good 
understanding of these networks. 

The properties of the small scale structure have been analyzed over the last 
decade and a half and its existence has deep implications for the networks and 
consequently for the cosmology. Strings with a fair amount of short distance struc- 
ture will self-intersect and produce small loops with high probability. However, 
the typical size at which loops form has been an issue of much debate over the 
years. Estimates range from the thickness of the string ^H] at the lower end, to 
sizes just an order of magnitude below the horizon distance Less radically, 

Refs. [381 EZj suggest that loops form predominantly with sizes 10^'^ — 10"^ times 

■^A notable exception occurs in superstring networks for which there are several attractors; 
which one is chosen at the outcome depends on the initial conditions in this case [40j . 
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smaller that the horizon. It should be noted that larger loops lead to enhanced 
signatures of several types, and consequently to tighter bounds on Gfi/c^ |47] . 

Of course, after establishing the presence of small scale structure, one would 
like to understand its evolution. Does the short distance structure approach a 
scaling regime as well? Part of literature assumes that gravitational radiation 
is needed for scaling, and that it determines the size of loops. However, the 
most recent simulations suggest that the late time properties evolve in a scaling 
fashion [IS], HSl ESI HB] , even though gravitational radiation is not included in 
their algorithms. 

1.4.3 Analytic approaches 

The present dissertation is based on analytic approaches to the study of scaling 
networks of cosmic strings. Several other studies have been previously conducted 
and we now briefly survey them. This will serve as a platform for comparison 
with the work presented in this dissertation. 

The one-scale model [35| B9| [50] nicely accommodates the scaling behavior of 
the network at late times as seen in simulations but only concerns the large scale 
properties. A model with two scales was considered in [5l], one of which char- 
acterized the long-string density and the other one being the persistence length. 
However, none of these explained the small scale structure observed in the sim- 
ulations and so a third length scale was later added |52j to specifically describe 
it. 

The picture arising from the detailed studies of [52j confirms that the problem 
of cosmic string network evolution is very complex: all the processes that occur 
during the evolution interact with each other. Unfortunately, the large number of 
unknown parameters (and assumptions) somewhat reduce the appeal of this three- 
scale model. The complexity of these systems was also revealed in attempts to use 
path- integral methods [521 El] • There, even working in fiat space and assuming 
the simplest possible string probability distribution computations quickly become 
quite involved. 

An improvement over the original one-scale model came with the velocity- 
dependent one-scale model [551 ES] , which allowed a varying averaged string ve- 
locity. This approach can address the evolution in transient regimes and therefore 
in realistic cosmological scenarios but, once again, ignores the small-scale issues. 

Along different lines, another model with a single scale was developed in |57] . 
The scale considered therein, determined by the number density of kinks, was 
found to approach very small values compared to the horizon distance but the 
model is in many ways too simplistic. 
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1.5 Outline 

The materials presented in the following chapters relie on Refs. |58l [591 EQl ED 
|62] and are mainly results of collaborations with Joseph Polchinski and Florian 
Dubath. 

As we have alluded to before, the evolution of cosmic string networks is a noto- 
riously difficult problem. Nevertheless, the problem is well-posed |63] and so our 
difficulties in studying such networks reflect their complexity as a whole and not a 
lack of understanding of the microscopies per se. Therefore, it is not unreasonable 
to imagine that a good strategy to tackle the evolution of such networks starts by 
focusing on a microscopic description. In Chapter 2 we set up an analytic model 
whose purpose is to describe the small scale structure on strings. The model is 
derived from a few simplifying assumptions on the network evolution which are 
expected to hold over a large range of length scales. This approach has the advan- 
tage of avoiding the nonlinearities of the system which dominate at the horizon 
scale. The full non-linear string equations do not seem amenable to analytic tech- 
niques and at such large scales we must trust the numerical simulations. The 
main conclusion is that the strings become smooth at short distances but there 
is a power law in this approach to straightness. The exponent x controlling the 
deviation from straightness is determined by cosmological parameters and large 
scale properties of the networks. 

In Chapter 3 we use our model to investigate the issue of gravitational radiation 
from long cosmic strings. The emission of gravitational radiation depends crucially 
on the spectrum of perturbations along the string and also leads to the smoothing 
of cosmic strings below a certain scale. We find that this scale is smaller than the 
horizon distance by a power of G^/c^ and the exponent is related to the spectrum 
of perturbations, namely x- 

Chapter 4 is devoted to a study of loop formation from self-intersections of 
cosmic strings and takes the model of Chapter 2 as input. The main surprise 
here is that without inserting an ultra-violet (UV) cutoff on the spectrum of 
fluctuations the rate of loop production diverges, despite the smoothness of the 
strings on small scales. A careful treatment reveals that loops of all sizes form 
simultaneously (instead of cascading down from large loops to small ones) with a 
power law distribution determined by x- This result is in good agreement with 
recent simulations [33, HO] . The role of the UV cutoff is taken by the gravitational 
radiation scale and it partly determines the normalization of the loop distribution. 

We reserve Chapter 5 to a survey of the small loop population properties and 
in particular to the investigation of the appearance of a scaling regime in the 
loop distribution. We point out that loops with sizes close to the gravitational 
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radiation cutoff are formed with large Lorentz boosts. Also, we show that the 
loop length distribution approaches a scaling regime at late times, both with and 
without the inclusion gravitational radiation. 

In Chapter 6 we consider some consequences for observational signatures. In a 
first part we calculate deviations from perfect lensing by cosmic strings and derive 
expectations for the alignment between several pairs of cosmic string lenses. Sub- 
sequently, we investigate the detectability of (quasi-) periodic gravitational waves 
generated by loops, taking into account the large Lorentz factors characteristic of 
small loops. We find that such mechanism of gravitational wave production can 
yield detectable signals for values of Gji/c^ of order (9(10~^), but this depends 
strongly on the direction of motion of the small loops relative to the detectors. 
When integrated over the population of loops within the horizon, this gives an 
expected rate of detectability at Advanced LIGO equal to 10~^ events per year, 
for Gn/c^ - 10-9. 

Finally, we conclude in Chapter 7 with some final remarks. 

In the remaining part of this dissertation we shall use units in which c—1. 
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Chapter 2 

Model for small scale structure 
on cosmic strings 

In this chapter we dehneate a model designed to describe the small scale struc- 
ture on cosmic strings. The basic object we consider is the two-point function 
which describes how points along a single string are correlated with each other. 
The full microscopic equations for the string network [531 ES] appear to be too 
complicated to solve, and so we model what we hope are the essential physical 
processes. In order to make the problem more tractable we resort to simplifying 
assumptions and well-motivated approximations. 

We begin by offering a preliminary Section [2TT] where the necessary ingredients 
from the cosmic string literature are introduced. In Section 12.21 we identify our 
assumptions, arguing that over a comprehensive range of scales the dominant effect 
in the evolution of cosmic string networks is the stretching due to the expansion of 
the universe. Based on these assumptions, in Section [231 we are able to determine 
the form of the two-point function up to two parameters which are inferred by 
matching our solution to numerical simulations. We will see that the string is 
actually rather smooth, in agreement with simulations [37j: its fractal dimension 
approaches one as we go to smaller scales. There is, however, a nontrivial power 
law that reveals itself in the approach of the fractal dimension to one. The critical 
exponent, which is determined by the mean string velocity, is related to the power 
spectrum of perturbations on the long string. Section 12^ is devoted to discussion 
and comparison of our results for the two-point function with the above-mentioned 
simulations. 

The model presented here was developed together with Joseph Polchinski in 
Ref. [58j and follows that reference closely. 
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2.1 Preliminaries: Basic ingredients 

In Chapter 1 we saw that cosmic strings can arise both in field theory (namely 
in Grand Unification Theories) and string theory contexts. Cosmic superstrings 
are honest one-dimensional objects but their field theoretic counterparts can also 
be regarded as one-dimensional defects, at least on scales much larger than their 
thickness. Thus, it comes as no surprise that for the simplest class of "vanilla" 
cosmic strings their dynamics is described by the Nambu action: 



A string evolved in time sweeps out a two-dimensional surface, the worldsheet. 
The above action is simply proportional to the area swept by the string in space- 
time. The time-like and space-like coordinates parameterizing the worldsheet are 
denoted by cr° and a^, respectively, and the integral is taken over these coordinates. 
The prefactor /i is nothing but the string mass per unit length and the quantity 
7 is the determinant of the metric induced on the worldsheet from the embedding 
in four- dimensional spacetime. Defining this embedding by the function^ a;'^ (a"), 
the induced metric is given by 



Expression ( 12. ip represents an effective action for gauge strings which can be de- 
rived under rather general considerations and provides a good description when- 
ever their radius of curvature is much larger than the thickness of the string [6]. 

The action (12. ip enjoys the property of invariance under a large set of trans- 
formations, namely reparametrizations of the worldsheet. This freedom can be 
used to bring the worldsheet metric into diagonal form. Furthermore, one can 
also identify the time-like coordinates of the worldsheet and of the ambient four- 
dimensional spacetime, 0"° = x° = t. Such a procedure still leaves a residual gauge 
freedom corresponding to time- independent reparametrizations of Therefore, 
in this so-called transverse gauge the spacetime coordinates of the string are fully 
specified by the spatial vector x(t, a), where a parametrizes the string at fixed time 
t. Note that the first gauge fixing condition amounts to the following constraint: 



"'^ Greek indices will be used for spacetime coordinates (/i — 0, . . . , 3) and lowercase Latin 
indices will represent worldsheet coordinates {a — 0,1). 

^ Since the time-like coordinate is already fixed in this gauge we shall drop the superscript 
from cr^ in the rest of this dissertation. 




(2.1) 



lab = 9f,udaX^dbX'' . 




(9tx ■ (9o-x = . 



(2.3) 
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This states that at any given point on the string the velocity is orthogonal to the 
tangent vector, thus justifying the name of this particular gauge. 

We are interested in cosmic strings evolving in cosmo logical spacetimes of the 
Friedmann-Robertson- Walker (FRW) type, for which the ambient metric Qf^i, is 
defined by 

^ g^^dx^dx" = -dt^ + a{tfdx ■ rfx = a{Tf{-dT'^ + rfx ■ rfx) . (2.4) 

The conformal time r has been introduced via dr = a~^{t)dt and the scale factor 
a{t) determines the expansion of these spatially flat universes. Therefore, two co- 
moving observers (sitting at fixed x coordinates) separated by a physical distance 
^ experience a recessional velocity relative to each other equal to f = HI, where 

a dt 

is the Hubble parameter, whose present value is about 73 Kms~^MPc~^. 

We shall consider only the situation of universes dominated by radiation or by 
matter, in which case the scale factor grows in time as a simple power law, 

a oc oc r'^' , u' = z//(l - u) , (2.6) 

where i/ = 1/2 and u' = 1 [u = 2/3 and u' = 2) in a radiation (matter) dominated 
era. This implies that the Hubble parameter decreases with time as t~^. There- 
fore, a useful measure for timescales is provided by the Hubble time, H~^. This 
represents the time it would take for the universe to expand to twice its (linear) 
initial size assuming the expansion rate remained fixed at the initial value. 

An important feature of FWR spacetimes is the presence of a cosmological 
horizon: events can only be causally connected if they are space-separated by no 
more than the horizon distance 

dH{t) = a{t) r^ = 7^. (2.7) 
Jo a{t') 1 - 

The second equality applies only to the case of a power law expansion. In as- 
trophysics the cosmological redshift z is commonly used to specify distances; it 
contains the same information as the scale factor since the two are related by 
1 -|- z{t) = a{t)~^, where conventionally the scale factor at present time is set to 
unity. 

In an FRW background the equation of motion governing the evolution of 
a cosmic string, derived by varying the action ( 12. ip with respect to the string 
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embedding, is [61] 




(2.8) 



Here e is given by 




) 



1/2 



e = 



(2.9) 



These equations hold in the transverse gauge, now with the time-hke worldsheet 
coordinate identified with the conformal time rjfl It is the non-hnear nature of 
these equations that undermine a complete analytic treatment. 
The evolution of the parameter e follows from equation (12.81) . 



and it is related to the energy dE carried by a segment of string with infinitesimal 
coordinate length da in an expanding universe by dE = fia{T)e{T,a)da. The 
denominator in equation (12. 9p neatly accounts for the Lorentz factor expected for 
relativistic motions. 

From the second derivative terms in equation (12.81) it follows that signals on 
the string propagate to the right and left with da = ±dT/e. The equation of 
motion also includes a friction term which is controlled by the Hubble parameter 
d/a. Thus, in a fiat spacetime the structure on a short piece of string at a given 
time is a superposition of left- and right-moving segments, and it is these that we 
follow in time in Section 12. 3[ In an expanding universe the left- and right-moving 
waves interact — they are not free as in fiat spacetime. 

Let us now turn to the issue of evolution of cosmic string networks. As we 
have discussed in the Introduction, a key property of such networks is that their 
energy density ps does not dominate over the radiation or matter energy densities 
in their respective epochs (for which pr oc ait)'"^ oc and pm oc a(t)^^ oc t^^). 
If the only process acting on the strings were the conformal expansion imposed 
by the growth of the universe we would have ps oc a(t)~2. However, there are 
several other mechanisms that come into play and reduce ps down to a (small) 
constant fraction of the background energy density when the scaling regime is 
reached in each era. Indeed, when scaling, the length of string within a horizon 
volume grows with t (so the energy density contained in the cosmic string network 
is proportional to pt~^) and the ratios Ps/Pr and Ps/Pm are proportional to the 

■^From now on we will use dots and primes to refer to derivatives relative to the conformal 
time T and the spatial parameter a along the string, respectively. 




(2.10) 
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small parameter Gj^. This fact is at the origin of the cosmological viability of 
cosmic strings. We now discuss the individual processes mentioned above. 

Stretching 

As the spacetime fabric expands, it tends to 'drag' with it all objects sitting 
within. As the FRW models are homogeneous and isotropic, one would naively 
expect that infinite strings would grow uniformly. This is true on scales larger 
than the horizon size dn for which irregularities on a string are just conformally 
amplified. However, considering small perturbations on an otherwise straight 
string and performing a linearized analysis [65] one finds that the modes with 
a wavelength shorter than dn do not grow in amplitude accordingly. Thus, on 
smaller scales the string effectively straightens. This effect is incorporated in the 
description of the string evolution by the Nambu action in curved spacetime. 

Gravitational radiation smoothing 

Another source of smoothing comes from gravitational radiation. Cosmic 
strings interact gravitationally through their energy-momentum tensor. Straight 
strings do not generate any gravitational effect on surrounding matter but this 
is no longer true once we introduce oscillations. In particular, any given point 
on the string interacts with the rest of the defect, resulting in power emitted as 
gravitational waves (GW). This in turn steals energy from the cosmic string itself 
which then consequently straightens. The strength of the GW emission is natu- 
rally controlled by the dimensionless string coupling G/i, the product of Newton's 
constant and the string mass per unit length. However, this smoothing by gravi- 
tational radiation becomes relevant only below a length scale proportional to dn 
and to a positive power of G/x P, [66l [671 [68]. We will have more to say about this 
in Chapter 3. 

Intercommuting 

As the cosmic string network evolves in time, often two segments of string 
intersect. When this happens and there is just one type of string in the network 
only two outcomes are possible: either the segments reconnect (with probability 
P) or they pass through each other (with probability 1 — P)Q This process 
is shown in Fig. 12. 1[ Field theory strings always reconnect, except when very 

'*These two classical processes dominate over situations resulting in entangled strings [69j . 
Also, In more complicated networks one might have Y-junctions as well and therefore 'bridges' 
between strings [70] , 
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Figure 2.1: The top panel shows an intercommutation event in which two segments 
of string exchange partners. The bottom panel displays a single cosmic string self- 
intersecting and chopping off a loop. In such processes pairs of kinks are formed 
on the intervening segments and travel in opposite directions along the string. 

extreme conditions are met by the collision parameters [69], whereas for cosmic 
superstrings the reconnection is a quantum process whose strength depends on 
the value of the string coupling constant [71j. 

An intercommutation event can result in the formation of closed loops. This 
important process can happen when a string self-intersects, as illustrated in the 
bottom panel of Fig. 12. 1[ A side effect of a reconnection is the appearance of 
pairs of kinks on each of the two intervening segments. These kinks then start 
moving apart and open up (becoming less acute) as time elapses, again due to the 
expansion of the universe. 

Decay of loops by emission of gravitational radiation 

Loops with sizes smaller than the horizon scale do not grow with the expan- 
sion of the universe. These loops live in an approximately flat spacetime and so 
they oscillate almost freely, instead. Consequently, they emit gravitational radia- 
tion and eventually decay completely, disappearing from the network. Thus, the 
formation of small loops and their subsequent decay are essential to achieve the 
all-important scaling regime. 
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2.2 Assumptions of the model 

We consider "vanilla" cosmic strings, a single species of local string without 
superconducting or other extra internal degrees of freedom. In this dissertation 
we will assume for simplicity that every collision between two segments of string 
results in a reconnection, so that the intercommutation probability is P = 1. The 
evolution of a network of such strings is dictated by three distinct processes!^ 

• First, the expansion of the universe stretches the strings; 

• Gravitational radiation also has the effect of straightening the strings, but it 
is significant only below a length scale proportional to dn and to a positive 
power of Gfi. Since this is parametrically small at small G^, we will ignore 
this effect for the present purposes; 

• Finally, intercommutations play an important role in reaching a scaling so- 
lution, in particular through the formation of closed loops of string. At first 
we shall neglect this effect but we will be forced to return to this issue in 
Chapter 4. 

Let us consider the evolution of a small segment on a long string. We take the 
segment to be very short compared to the horizon scale, but long compared to the 
scale at which gravitational radiation is relevant. The scaling property of the net- 
work implies that the probability per Hubble time for this segment to be involved 
in a long string intercommutation event is proportional to its length divided by 
(Ih, and so for short segments the intercommutation rate per Hubble time will 
be small. Formation of a loop much larger than the segment might remove the 
entire segment from the long string, but this should have little correlation with 
the configuration of the segment itself, and so will not affect the probability dis- 
tribution for the ensemble of short segments. Formation of loops at the size of the 
segment and smaller could affect this distribution but this process takes place only 
in localized regions where the left- and right-moving tangents are approximately 



equalO Thus, there is a regime where stretching is the only relevant process. 

If we follow a segment forward in time, its length increases but certainly does 
so more slowly than the horizon scale dn, which is proportional to the FRW time 
t. Thus the length divided by dn decreases, and therefore so does the rate of 

^We neglect damping forces arising from interactions with the surrounding radiation back- 
ground. Such frictional terms enter the evolution equations in the same manner as the Hubble 
damping and become subdominant roughly after a time (G/i)~^if , where if is the time of string 
formation [6]. 

^Nevertheless, the results of Chapter 4 indicate that the production of small loops is large. 
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intercommutation. If we follow the segment backward in time, its length eventu- 
ally begins to approach the horizon scale, and the probability becomes large that 
we encounter an intercommutation event. Our strategy is therefore clear. For 
the highly nonlinear processes near the horizon scale we must trust simulations. 
At a somewhat lower scale we can read off the various correlators describing the 
behavior of the string, and then evolve them forward in time using the Nambu 
action until we reach the gravitational radiation scale. The small probability of 
an intercommutation involving the short segment can be added as a perturbation. 
This approach is in the spirit of the renormalization group, though with long and 
short distances reversed. 



2.3 Two-point functions at short distance 

From Eq. f l2.10p it follows that the time scale of variation of e is the Hubble 
time, and so to good approximation we can replace with the time-averaged v'^ 
(bars will always refer to root-mean-square (RMS) averages), giving e oc a~^^ as 
a function of time only0 The RMS velocities for points on long strings are taken 
from simulations 

radiation domination: vl ^ 
matter dommation: t> ^ 

From the definition of e it follows that the energy of a segment of string of 
coordinate length a is E = jia{T)e{T)a. For simphcity we will refer to E / n as the 
length / of a segment, 

/ = a{T)t{T)a , (2.12) 

though this is literally true only in the rest frame. 
As a consequence, 

/ocr^'ocf^, C' = (1 - 2t;^)i^' , C = (1 - 2t;^)i^ • (2.13) 

In the radiation era ~ 0.1 while in the matter era Cm ~ 0.2. Thus the phys- 
ical length of the segment grows in time, but more slowly than the comoving 
length ^65], and much more slowly than the horizon length dn- 

For illustration, consider a segment of length (10~^ to 10~^)(i//, as would 
be relevant for lensing at a separation of a few seconds and a redshift of the 

^The transverse gauge choice leaves a gauge freedom of time-independent a reparameteriza- 
tions. A convenient choice is to take e to be independent of a at the final time, and then e will 
be (j-independent to good approximation on any horizon length scale in the past. 
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order of z ~ 0.1. According to the discussion above, l/dn depends on time as 
t'^~^ ~ in the matter era. Thus the length of the segment would have been 
around a hundredth of the horizon scale at the radiation-to-matter transition. In 
other words, it is the nonlinear horizon scale dynamics in the radiation epoch 
that produces the short- distance structure that is relevant for lensing today, in 
this model. This makes clear the limitation of simulations by themselves for 
studying the small scale structure on strings, as they are restricted to much smaller 
dynamical ranges. 

For our purposes it is convenient to work not with the velocity and tangent 
vectors but instead with a linear combination of them, p-t = x ± -x'. In the 
transverse gauge we are adopting, these are unit vectors as can easily be seen 
by employing the definition (12. 9p . Furthermore, in a flat spacetime (d = 0) the 
vector p_|_ (p-) would be a left-mover (right-mover) as it is annihilated by the 
differential operator dr — -d^- In terms of these left- and right-moving unit vectors 
the equation of motion (12. 8p can be written as [10] 

p± T -P± = -- [Pt - (P+ ■ P-) P±] • (2.14) 

We will study the time evolution of the left-moving product p+(r, a) ■ p+(r, a'). 
For this it is useful to change variables from (r, a) to (r, s) where s is constant 
along the left-moving characteristics, s — s'/e = 0. Then 

5r(p+(s,r) ■ p+(s',r)) = --(p_(s,r) • p+(s',r) + p+(s,r) ■ p_(s',r) 

a \ 

+a{s, r) p+(s, r) ■ p+(s', r) + a{s', r) p+(s, r) ■ p+(s', r)^ , (2.15) 

where a = — p+ ■ p_ = 1 — 2f ^. 

The equations of motion (I2.141[2TT5|) are nonlinear and do not admit an analytic 
solution, but they simplify when we focus on the small scale structure. If p+(s, r) 
were a smooth function on the unit sphere, we would have 1 — p+(s, r) ■p_|_(s', r) = 
0{[s — s']^) as s' approaches s. We are interested in any structure that is less 
smooth than this, meaning that it goes to zero more slowly than [s — s']^. For 
this purpose we can drop any term of order [s — s']^ or higher in the equation of 
motion (smooth terms of order s — s' cancel because the function is even) . 

Consider the product p_(s, r) ■ p+(s',r). The right-moving characteristic 
through (s, r) and the left-moving characteristic through [s', r) meet at a point 
(s, T — S) where S is of order s — s'^ Eq. (I2.14p states that p+ is slowly varying 

^Explicitly, for given r we could choose coordinates where e(T) = 1 and s{t' ,a) — a + t' — 
T + 0([r' - r]2), and then S = {s - s')/2. 
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along left-moving characteristics (that is, the time scale of its variation is the FRW 
time t), and p_ is slowly varying along right-moving characteristics. Thus we can 
approximate their product at nearby points by the local product where the two 
geodesies intersect, 

p_(s, r) • p+(s', r) = -a(s', r-5) + 0(s- s') . (2.16) 

Then 

dr(p+(s,T) ■p+(s',t)) = -(a(s',r - 5) + a{s,T + 5) 

a \ 

-a{s, t) p+(s, t) • p+(s', r) - a{s', r) p+(s, r) • p+(s', r)^ + 0{s - s') . (2.17) 

When we integrate over a scale of order of the Hubble time, the S shifts in the 
arguments have a negligible effect 0{S) and so we ignore them. Defining 

h++{s, s', r) = 1 - p+(s, r) • p+(s', r) , (2.18) 

we have 

drh++{s, s', t) = --/i++(s, s', t)[q;(s', t) + a{s, r) + 0{s - s')] . (2.19) 
a 

Thus 

/i++(s,s',Ti) = /i++(s,s',To)exp|-i/'y" ^[a(s',r) a(s,r) 0(s - s')] • 



(2.20) 

Averaging over an ensemble of segments, and integrating over many Hubble 
times (and therefore a rather large number of correlation times) the fluctuations 
in the exponent average out and we can replace q;(s, r) with a — 1 — 2v^, 

(/i++(s,s',ri)) ^ (/i++(s,s',To))(ri/ro)-^'^'^ . (2.21) 

Note that in contrast to previous equations the approximation here is less con- 
trolled. We do not have a good means to estimate the error. It depends on the 
correlation between the small scale and large scale structure (the latter determines 
the distribution of a), and so would require an extension of our methods. We do 
expect that the error is numerically small; note that if we were to consider instead 
(ln/i++(s, s', r)) then the corresponding step would involve no approximation at 
all. 
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Averaging over a translationally invariant ensemble of solutions, we have 

(/i++(a - a', n)) ^ (ri/ro)-2^'^(/i++(a - a', Tq)) . (2.22) 

We have used the fact that to good approximation (once again in the sense of 
Eq. (12.2ip ). e is only a function of time, and so we can choose cr = s — J dr/e and 
a — a' = s — s'. Equivalent ly, 

{h+Aa-a',r))^l^^^. (2.23) 

The ratio of the segment length to dn = + is 

I aeia — a') a — a' . . 

— oc ^ oc . 2.24 

The logic of our earlier discussion is that we use simulations to determine the 
value of h^^ at l/dn somewhat less than one, and then evolve to smaller scales 
using the Nambu action. That is, 

h++{a - a', t) = ho when a - a' = Xot^+^''''"' , (2.25) 

for some constants xq and ho. We assume scaling behavior near horizon scale, 
so that Hq is independent of time. Using this as an initial condition for the 
solution (12.231) gives 

{h^4cr-a',r))^ho( ^~X,) \ ^ A{l/ty\ (2.26) 



where 



1 + lu'v^ \ — va 

In the last form of equation (12.261) we have expressed the correlator in terms of 
physical quantities, the segment length / defined earlier and the FRW time tj^ 
Eq. (I2.26P is our main result. Equivalently (and using o parity), 

(p+(a, r) ■ p+(a', r)) = (p_(a, r) ■ p_(a', r)) ^ 1 - Ailjtf^ . (2.28) 

In the radiation era Xr ~ 0.10 and in the matter era Xm ~ 0.25. The most 
relevant constants that determine the evolution of small scale structure on cosmic 



^We have not yet needed to specify numerical normalizations for a and t, or equivalently for 
e and a. The value of ho depends on this choice, but the value of A does not. 
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Cosmological epoch 


V 




a 


X 


Radiation dominated 


1/2 


0.41 


0.18 


0.10 


Matter dominated 


2/3 


0.35 


0.30 


0.25 



Table 2.1: Summary of relevant quantities for the evolution of small scale structure 
on cosmic strings in a scaling regime for both radiation and matter dominated eras. 
All quantities are defined in the body of the text. The RMS velocities were 
taken from simulations by C. Martins and E. P. S. Shellard (2005). 

strings in a scaling regime are summarized in Table 12.11 for both the radiation- 
and matter-dominated eras. 

There can be no short distance structure in the correlator p_,_ ■p_, because the 
left- and right-moving segments begin far separated, and the order aja interaction 
between them is too small to produce significant non-smooth correlation. Thus, 
from fl2.16l) we get 

(p+(or, r) ■ p_(a', r)) = -« + 0{g - g') . (2.29) 

2.3.1 Small fluctuation approximation 

Before interpreting these results, let us present the derivation in a slightly 
different way. The exponent x is positive, so for points close together the vectors 
p+((T, r) and p+((j',r) are nearly parallel. Thus we can write the structure on 
a small segment as a large term that is constant along the segment and a small 
fluctuation: 

p+(a, r) = P+(r) + w+(r, a) - ^P+(rX(r, a) + . . . , (2.30) 

with Pj^{tY = 1 and P+(t) ■ w_|_(r, cr) = 0. Inserting this into the equation of 
motion (12.141) and expanding in powers of w+ gives 

P+ = --[P--(P+-P-)P+] , (2.31) 

w+--w; = -[(w+.P_)P+ + (P+.P_)w+] 
c a 

= -(w+-P+)P+ + -(P+-P_)w+ . (2.32) 
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Since the right-moving p_ is essentially constant during the period when it crosses 
the short left-moving segment, we have replaced it in the first line of ( 12.32^ with a 
(T-independent P_. In the second line of (12.321) we have used w+ ■ P+ = 0. In the 
final equation for w+, the first term is simply a precession: P+ rotates around an 
axis perpendicular to both P_|_ and P_, and this term implies an equal rotation 
of w+ so as to keep w+ perpendicular to P+. Eq. (I2.32p then implies that in a 
coordinate system that rotates with P+, w+ is simply proportional to a~". 
It follows that 



scales as a~^" as found above (again we are approximating as in eq. (12.211) . and 
again this statement would be exact if we instead took the average of the log- 
arithm). Similarly the four-point function of w_|_ scales as a~^". We have not 
assumed that the field w+ is Gaussian; the n-point functions, just like the two- 
point function, can be matched to simulations near the horizon scale. We can 
anticipate some degree of non-Gaussianity due to the kinked structure; this will 
be discussed further in Chapter 6. 

2.4 Discussion 

Now let us discuss our results for the two-point functions. We can also write 
them as 



These are determined up to two parameters and A that must be obtained 
from simulations. A first observation is that these expressions scale as they are 
functions only of the ratio of / to the horizon scale. This is simply a consequence 
of our assumptions that the horizon scale structure scales and that stretching is 
the only relevant effect at shorter scales. We emphasize that these results are for 
segments on long strings; we will discuss loops in Chapter 4. 

It is natural to characterize the distribution of long strings in terms of a fractal 
dimension. The mean squared spatial distance between two points (also known 




(2.33) 




(2.34) 
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as the extension) separated by coordinate distance a is 

{r\l, t)) = (x' ■ x') f dl' f dl" corr,(/' - t) 

Jo Jo 



{2x+mx + 2){l-v 



(2.35) 



We can then define the fractal dimension df (which is 1 for a straight fine and 2 
for a random walk), 

d, - ^^^"^ ^ 1 + -^^(^A)'' + 0((l/t)'-) (2 36) 

- dHr^il, t)) ^ ^ + (2x + l)(2x + 2)(1 - + ^^^^^'^ ^ ' ^^.36) 

The fractal dimension approaches 1 at small scales: the strings are rather smooth. 
There is a nontrivial scaling property, not in the fractal dimension but rather in 
the deviation of the string from straightness, 

1 - corr^.(/, t) oc {l/tf^ , 1 - corrt(/, t) oc {l/tf^ . (2.37) 

We define the scaling dimension dg = 2%. Note that dg is not large, roughly 0.2 
in the radiation era and 0.5 in the matter era, so the approach to smoothness is 
rather slow. 

One might also consider the fractal dimension of the p± curves which live on 
the unit sphere. To this end we consider the extension between two such vectors 
separated by an infinitesimal worldsheet coordinate 6a, i.e. we take p+(cr) and 
p+(cr + 6a) ~ p+icr) + 5p+(cr). Now using our result fl2.28p we see that 

(p'+(^) ■ pV(^ + « • (2-38) 

Therefore, the extension is given by 

{\6p+\) = {M6a ^ 6a^ . (2.39) 

It then follows that the fractal dimension of the p± curves is l/x- Note that this 
exponent is quite large, in contrast with the smoothness of the string configuration 
itself; it equals 10 in the radiation dominated era and 4 in the matter dominated 
era. 

Our general conclusions are in agreement with the simulations of Ref. [37], 
in that the fractal dimension approaches 1 at short distance. To make a more 
detailed comparison it is useful to consider a log-log plot of 1 — corr2;(/,t) versus 
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10' 



X = 0.10 , A = 0.60 




10-^1 . . . . . . . I 

-3 -2 -1 1 2 3 4 5 

log(//^) 

Figure 2.2: Comparison of the model (dashed hne) with the data provided by 
C. Martins and E. P. S. Shellard (sohd red hne) in the radiation-dominated era, 
for which the correlation length is ^ ~ 0.30t. 

/, as suggested by the scaling behavior (12.371) F"! The comparison is interesting. 
At scales larger than dn the correlation goes to zero0 Rather abruptly below 
dn the slope changes and agrees reasonably well with our result. It is surprising 
to find agreement at such long scales where our approximations do not seem very 
precise. On the other hand, at shorter scales where our result should become 
more accurate, the model and the simulations diverge; this is especially clear at 
the shortest scales in the radiation-dominated era (Fig. 12.21) . Note that at these 
smaller scales the simulations seem to indicate a larger exponent x ~ 0-5, which 
corresponds to the functions p± mapping out a random walk on the unit sphere. 

One possible explanation for the discrepancy would be transient behavior in 
the simulations. We have argued that the structure on the string is formed at the 
horizon scale and 'propagates' to smaller scales (in horizon units) as the universe 
expands. In Ref. [37] the horizon size increases by a factor of order 3, and so even 
if the horizon-scale structure forms essentially at once, the maximum length scale 
over which it can have propagated is 3^"'', less than half an order of magnitude. 
At smaller scales, the small scale structure seen numerically would be almost 
entirely determined by the initial conditions. On the other hand, the authors 

^°We thank C. Martins for replotting the results of Ref. [37^ m this form. 
^^In terms of the correlation length ^ which is used in reference [37] the horizon scale is ~ 6.7^ 
in the radiation era and ~ 4.3^ in the matter era. 
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X = 0.25 , A = 0.61 
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Figure 2.3: Comparison of the model (dashed hne) with the data provided by 
C. Martins and E. P. S. Shellard (sohd blue line) in the matter-dominated era, for 
which the correlation length is ~ 0.69t. 

of Ref. [33| (private communication) argue that their result appears to be an 
attractor, independent of the initial conditions, and that loop production may 
be the dominant effect. Motivated by this we have examined loop production 
and these studies are presented in Chapter 4. Indeed, we find that this is in 
some ways a large perturbation, and it is conceivable that the chopping-off of 
small loops caused by self-intersecting long strings could significantly smoothen 
the latter. However, it still remains an open question whether this mechanism 
can resolve the discrepancy seen in the short distance structure of the two-point 
function. 

Thus far we have discussed corr^^. Our result fl2.34p implies a linear relation 
between corr^; and corr^. In fact, this holds more generally from the argument 
that there is no short- distance correlation between p+ and p_, Eq. (12.291) : 

(1 — 'i;^)corra;(/, t) — v'^corit^l, t) = 



Inspection of Fig. 2 of Ref. ^37j indicates that this relation holds rather well at all 
scales below ^. 

The small scale structure on strings is sometimes parameterized in terms of 
an effective mass per unit length [721 [73] . The basic idea is to consider a coarse- 




(2.40) 
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grained description of a cosmic string at some scale. Any wiggles on it with 
wavelengths much smaller than this scale will appear as smooth at the cost of 
introducing an effective mass per unit length /Xeg (and an effective tension as 
well). For a segment of length / the effective mass per unit length is given by 



where we have made use of result fl2.35p . Note that this is strongly dependent on 
the scale / of the coarse-graining. 

In conclusion, let us emphasize the usefulness of the log- log plot of 1 — corr^,. 
In a plot of the fractal dimension, all the curves would approach one at short 
distance, though at slightly different rates. The difference is much more evident 
in Figs. 12.21 and 12. 3[ and gives a clear indication either of transient effects or of 
some physics omitted from the model. 




(2.41) 
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Chapter 3 

Gravitational radiation from 
cosmic strings 

The "vanilla" cosmic strings we consider interact only gravitationally. A 
straight cosmic string leaves the surrounding space locally flat, merely introducing 
a deficit angle (giving rise to a conical space, globally) [25j. This is of course not 
the case once we allow for oscillating strings. In particular, gravitational radiation 
is emitted when waves traveling in opposite directions along the string collide and 
this effect is proportional to the dimensionless string tension Gfi. The resulting 
radiation depends on the power spectrum of the perturbations on the string [68] 
and so we can use the results obtained in Chapter 2 to address this matter. This 
is precisely what we present in this chapter. 

We begin by reviewing in Section 13.11 the computation of the power emitted 
in gravitational radiation from colliding small perturbations on long strings and 
the traditional picture of the gravitational back-reaction cutoff on the scales of 
fluctuations. In Section 13.21 this calculation is improved by incorporating the 
fact that colliding modes on a string only radiate efficiently if their wavelengths 
have the same magnitude and by employing the spectrum of perturbations found 
in Chapter 2. This gravitational back-reaction introduces a new scale in the 
network, below which the strings become smoother. Therefore, the two-point 
function obtained in the previous chapter must be corrected at small scales and 
this is done in Section 13. 3[ For later reference, we briefiy discuss the interplay 
between small scale structure and cusp formation in Section 13.41 

This chapter is mainly based on Ref. [59], with Joseph Polchinski. 
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3.1 Radiation from long strings and back-reac- 
tion 

As we mentioned in the introductory chapter, it is generally believed that, 
after their formation, cosmic string networks evolve into a scaling regime where 
all length scales grow linearly with the Hubble time. This would imply that 
the typical size of cosmic string loops is a fixed fraction a of the horizon scaleQ 
However, there is no consensus on the value of the proportionality constant a; in 
fact, estimates range over tens of orders of magnitude. 

A key question is whether the straightening of the string by gravitational 
radiation is necessary for scaling of the loop sizes. If so, a will depend on the 
dimensionless string tension Gfi. If not, then the evolution of the network will 
be purely geometric, and a could be a pure number, independent of Gfi. In 
Chapter 5 we will demonstrate that gravitational radiation effects are not needed 
to obtain scaling of the loop length distribution. On the other hand, in Chapter 4 
we will find that employing the simple stretching model down to arbitrary scales 
results in a UV divergence in the loop production function. Since smoothing by 
gravitational radiation is after all a real effect, it actually determines the cutoff 
and so we argue that a does in fact depend on Gfi. 

For many years, simulations appeared to show loops forming at the short 
distance cutoff scale, and this was interpreted as implying the need to include 
gravitational radiation. The necessity to consider it in order to achieve scaling 
of small scale structure was also suggested by the detailed study of Ref. [52] . 
However, the conclusions of [571 ES] were opposite: those authors found that 
stretching by the expansion of the Universe alone is sufficient to guaranty scaling 
of the small scale structure. Indeed, several recent simulations, indicate that the 
loops are actually forming above the cutoff scale. Refs. [15], |16] find notably large 
loops, a ~ 0.1, but this is superposed on a power law distribution that grows 
toward smaller scales. Refs. [Ml EZj find loops a few orders of magnitude smaller, 
but the distributions are still evolving in a non-scaling fashion. None of these 
simulations include gravitational radiation directly. 

When a perturbation on a string encounters another one traveling in the op- 
posite direction, gravitational waves are generated. The energy radiated away is 
stolen from the colliding fluctuations and as a result the string straightens. How- 
ever, this smoothing of the long strings by gravitational radiation is the dominant 
effect only below some critical length scale, usually called the gravitational radi- 
ation scale. This (time-dependent) scale then acts as a lower cutoff on the sizes 

^ There is not complete agreement even on this point: see Ref. j74j and references therein. 
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of loops produced at any given time. It was long assumed that gravitational ra- 
diation smoothed the long strings on scales of order the horizon length times Gfi, 
so this would be the size of loops if they formed at the gravitational radiation 
cutoff [751167]. However, gravitational radiation is an efficient energy loss channel 
only when the perturbations have comparable wavelengths [76]. Indeed, those 
authors showed that the radiation from long strings had been overestimated, and 
so it becomes important only at a shorter length scale, proportional to a larger 
power of Gfi [BB] ; the exponent depends on the power spectrum of fluctuations on 
the long string. 

The model developed in Chapter 2, which determined the form (12.281) for the 
two-point function at short distances in a scaling regime, should then provide 
some answer for the gravitational radiation scale. In the small fluctuation ap- 
proximation, the short distance structure can be equivalently expressed in terms 
of the fluctuations w-t as 

l-{p±{^,t)-p±{a',t)) = i([w±(a,t) - w±(a',t)]^) +0«) ^ Ail/t)'^ . (3.1) 

The model also predicts that the exponent x which controls the deviation of the 
string from straightness is determined by the RMS velocity v'^ of the long string 
population. We have noted that at small scales the simulations of [37] disagree 
from our model and seem to approach a larger value for the exponent, namely 
X ~ 0.5, corresponding to the functions p-t mapping out a random walk on the 
unit sphere. In any event, our discussion below just uses the general power law 
form (12.281) . and one can insert any assumed values for x and A. 

Gravitational radiation is important on scales short compared to the Hubble 
time and so we can use the flat metric g^j^i, = diag(— 1, +1, +1, +1), and also make 
the gauge choice ae = lH The modes are then functions only of -u = t + cr or 
V = t — cr. This is not surprising, given that the string equations of motion reduce 
to a simple wave equation in the coordinates t and a for a flat spacetime. The 
solution is a superposition of left- and a right-moving waves, 

x^'{u,v) = -[a''{u) + b''{v)] . (3.2) 

Hence, we may write 

(p-(^)-P-K)) ^ l-A[iv-v')/t]'^- (3.3) 

^This choice is possible when we consider scales small compared to the Hubble time. In 
Section [3.31 when we again consider cosmological evolution, we must reintroduce a{t) and e(t). 
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(the explicit t on the right-hand side is effectively a fixed parameter, varying only 
over the Hubble time). 

The power radiated in gravitational waves by an infinitely long source was 
first considered in ^77j. By employing a weak- field approximation, the gravita- 
tional field can be easily related to the energy-momentum tensor of the string and 
then the energy flux through a large radius cylinder surrounding it gives the power 
radiated, proportional to Gfi'^. A general formula for the gravitational radiation 
was obtained in [BB] under the assumption that the colliding waves traveling along 
the string have small amplitudes. Essentially, gravitational radiation is only emit- 
ted when left- and right-moving perturbations interact. This is manifest in the 
expression for the energy radiated per unit solid angle in the k direction, given 
below in the notation of 1761, 

dAE G^i^ 



POO 

/ du iu'' {\A\^\B\^ + \A* ■ B\^ - \A ■ B\^} , (3.4) 
Jo 



dVt 167r2 

with the left- and right-moving contributions given respectively by 
A^{k) = I dup^{u) exp <^ —- j du k ■ p+{u] 



oo 

+ 00 



B^ik) = j dvp''_{v)exp}^-^ j dv'k-p4v')^ . (3.5) 

Here, k^^ is the gravitational radiation wave vector, and p^. are null 4- vectors whose 
time component is identically 1 and with p± being the unit vectors aboveEl In 
the spirit of the small fluctuation approximation we orient a given segment of 
string mainly along the z-axis and consider fluctuations w-|- in the perpendicular 
xy-plsnae. Thus, 



P+{u) ^ (l,w+{u),^l-wl{u) 



p'^.iv) ~ \^l,w^{v),-^l-wUv)j . (3.6) 

We can ignore the last two terms in equation (13. 4p : they are equal when B'^ is 
real, and so cancel when ensemble-averaged over the short distance structure. 

Let us first review the argument of Refs. [76l[68]. First linearize the modes (13.51) 
in the oscillations, so the exponential factors become e*'^"" and e*'^''^ respectively, 
with 

uj = ka + kb , kz = kb - ka . (3.7) 



^These are related to the standard notation via p'^(u) ~ a^^'{u) and p^{v) = h^^'{v). 
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Here ui and kz are the frequency and wavenumber of the radiation. Then the 
ensemble averages in the hnearized approximation (^) are 



/ + 00 
dne*^""-^l"l ((p+(n) • p+(0)) - 1) = LAc^k-^'+'^k^'^ 
-oo 

'(l+2x).-2x 



{\B{k)\') ^ hAc^kp'^'''>t-'^ . (3.8) 

Here = 2 sin(7rx) r(l + 2%), taking the values 1.25 and 0.57 in the matter- and 
radiation-dominated eras respectively. The convergence factor e""^'"' accounts for 
the decay of the correlations on horizon scales. Its detailed form is unimportant: 
it only affects the Fourier transform for ka of order the inverse horizon size, while 
for gravitational radiation the relevant ka is much larger. The factors and ly are 
volume regulators: we artificially cut the oscillations off to produce finite trains, 
but these regulators of course drop out when we consider rates per unit time 
and length. Using dudcosO = 2dkadkb/{ka + kb), we can write the total energy 
radiated as 

(^,,2 roo POO 

{AE) = dkj dhika + h){\A{ka)f){\B{h)f) 

47r Jo Jo 

- ^^'^'''H.X Hdha Hdh „ ■ (3.9) 



4vr ^'"^Vo "Jo "{kahy+'H^^ 

The volume of the world-sheet is = ^lu^v, so this translates into a power per 
unit length 

dh-——r-r—. (3.10) 



dz/ 271 7o Jo (kahy+^H^x 

The total energy of the wavetrains is [66] 

poo ji, poo jr, 

E^^/ ^|A(Up + /./ ^\Bikb)\'. (3.11) 

Jo Jo 

Isolating a momentum range dka, we have 

A\Aika)\' ^ ~^\A{ka)\' r dh{\B{h)\^) . (3.12) 
^ Jo 

We have used the fact that momentum conservation determines that the energy 
coming from the A and B modes is in the proportion ka to k^. A given point u 
interacts with the v train for a time \ly, so the rate of decay becomes 

d r°° dk 

j^\A{ka)\' ^ -Gf,ka\A{ka)\'Ac^ . (3.13) 
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The integral is dominated by the lower hmit: most of the energy loss from the 
left-moving mode ka comes from its interaction with right-moving modes of much 
longer wavelength. Cutting off the lower end of the kb integral at the horizon 
scale, fcfe ~ we find that the decay rate is of order Gfika- This is faster than 
the Hubble time for 

ka > 0{GiJ,t)-^ , (3.14) 

meaning that modes with wavelengths A < OiGjJit) are exponentially suppressed. 
Thus, Eq. f l3.14p reproduces the usual (naive) estimate for the gravitational damp- 
ing length, namely that fluctuations on scales smaller than ~ G/it quickly decay 
due to gravitational radiation. 

3.2 Radiation and back-reaction improved 

We have already mentioned that including the fact that a perturbation on a 
string does not interact with the same efficiency with all other modes leads to a 
shorter gravitational radiation length scale [TBI [68] . Here we rederive these results 
in perhaps a simpler way. Our result for the exponent of G/i differs somewhat from 
Ref. [68]. This difference arises in going from the periodic wave train calculation 
of Ref. [76] to the random distribution on the actual network. Also, we use input 
from the model developed in Chapter 2 and from simulations [37] to determine 
the actual power spectrum on the long string. 

The calculation we carried out in the previous section implies that the energy 
loss from the left-moving modes ka comes primarily from their interaction with 
much longer right-moving modes at the horizon scale. Ref. [76] argued that this 
was paradoxical, because a wave that encounters an oncoming perturbation with 
much longer wavelength is essentially traveling on a straight string, and such a 
wave will not radiate [78] . They showed that this paradox arose due to the neglect 
of the exponential terms in the modes (13.51) . With their inclusion, the radiation 
is exponentially suppressed when the ratio of ka to kh becomes sufficiently large. 

Ref. [76] considered monochromatic wavetrains 

w+(m) = taCos{kau) , w''_{v) = ebCos{kbv) , (3.15) 

and found that keeping the full form (13. 5p of the amplitude gives a large suppres- 
sion when 

h/el<ka or kb>ka/el. (3.16) 

The first of these relations will cut off the integral (I3.13P below the horizon scale, 
but first we need to extend it to the incoherent spectrum on a long cosmic string. 
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We will do this in a systematic way below, but we can anticipate the answer. For a 
continuous spectrum, of course, a single frequency makes a contribution of measure 
zero. Indeed, the units are wrong: the Fourier transform of the fluctuation w_(?;) 
obeys 

(w_(A;) ■ w^k'Y) = 2'K5{k - k')e^{k) , (3.17) 

where e^{k) oc k~^^~^'^^H~'^^ has units of length, whereas in eq. (I3.16P the ampli- 
tudes must be dimensionless. Thus, in going to the continuum, we must replace 

el ^ k,e\k,) ^ {ht)~^^ . (3.18) 

The amplitude is then suppressed unless 

(M)i+2x>A;,t . (3.19) 

Using this lower cutoff in the decay rate ( 13.13^ gives 

^\A{ka)\' ~ -G^^{katf/^'^'^h-'\A{ka)\' = 

at 

This is faster than the Hubble rate for 

A;,>(G^)-(i+2x)^-i _ (3 21) 

The gravitational length scale is reduced from the naive 0{Gfit) by a factor 
{G^y^. The qualitative conclusion is the same as in Ref. [68], but the expo- 
nents do not seem to agree. The result there was ka > {Gfx)~^^^'^'^^^'^t~^ . It is not 
clear whether /? is to be identified with x + 1 (^^s suggested by Eq. (21) of Ref. |68j, 
compared with Eq. (I3.13P above) or with x (as suggested by [68] Eq. (31)), but in 
either case the exponents differ. This stems from a different method of convert- 
ing from the single-mode result to the continuous spectrum; our exponent will be 
borne out by the more formal treatment below. 

Finally, our analytic model of Chapter 2 gives for the exponent 1+2% the values 
1.2 in the radiation era and 1.5 in the matter era; using instead the simulations [37] 
would yield an exponent around 2.0 at the shortest scales in both eras0 

Numerically, there are two changes from Ref. |68j : the corrected expression for the exponent, 
and a more accurate estimate of the power spectrum. Ref. |68j effectively estimated the latter 
using = 0, so that their Eq. (31) is (3 = x\v^=o — ^/(l ~ v)- 



(3.20) 
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The necessary correction to the naive radiation formula comes from the previ- 
ously neglected exponential factor in B^{k), Eq. (13. 5p [76]. Expanding the expo- 
nential to second order in the fluctuations gives 



dv / dv' ip4v) ■ p4v') - 1] 

-oo J — oo 



exp i^tktiv' -v)-'^- j"^ dv"w4v") dv"wl{v")^ . (3.22) 

Noting that k\ = 4:kakb, when the first term in the exponent is of order one, the 
second and third terms in the exponent are respectively of order (katy^"^ {kbt)~^^'^~^ 
and kat{kht)~^~'^^. Thus, as k^ decreases they become important at precisely the 
scale (13.191) where the small-amplitude calculation is expected to break downjf] 
The linearized form for A{k) remains valid, so we can simply insert the corrected 
form (13:221) into the decay rate (l312|) . 

By completing the square, the exponent in Eq. ( 13.221) can be written as 

/ , \ 2 

iv' -V) I ka 



4:ka \ \/v' - V Jy 

ikz 



dv"wjv") 



8{v' - v) 



r dv" r dv'" [w^v") - w_(t;'")]' . (3.23) 

J V J V 



Converting j dkb —>■ {Aka) ^ J dk^ —>■ (ATcka) ^ J d'^k±, we have at fixed ka the 
integral 



poo p-\-oo p-\~oo 

/ dh\B{k)\'^ = dv dv' 

Jo J —CXI J —oo 

'^^{-W^) /" 'd«'"|w_K) -w_(,/")l^} .(3.24) 

We have used —k^ ^ ka ^ kf,, since this is the regime where the correction is 
important. 



^The terms compete at this scale because the smaU fluctuation w_ {k'^^) and the small number 
{kb/ka)^^'^ become comparable. Further terms in the exponent, which have been dropped, are 
of higher order in the small parameter. This is in keeping with Ref. [76], which also found that 
only terms up to quadratic order in the exponent were relevant. 
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We now approximate, replacing both the exponent and the prefactor with their 
mean values from equation ( 13. ip . Then 



dh{\B{k)\'^) = 2AlAm 



dv 



exp 



ikaA 

JO 

-^2x/(l+2x) 



dv" / dv 



Here 



2 sin 



2X 



r 



2X 



[4(1 + 



2x/{l+2x) 



1 l/(l+2x) 



l + 2x. 



(3.25) 



(3.26) 



is 2.5 in the matter era and 2.1 in the radiation era. We can improve this ap- 
proximation if we make the assumption that the ensemble is Gaussian. The first 
order correction, keeping contractions between the prefactor and the exponent, in- 
volves a straightforward calculation; its inclusion simply renormalizes the constant 



^1 + ex)&x' with 



x(i + x) 



r(i + 2x)^ 1 

[r(2 + 4x) 2(l + 4x)J ■ 



(3.27) 



However, this represents a small correction since = 0.035 in the radiation epoch 
and = 0.045 in the matter epoch. Thus, to good approximation. 



1 



b^Gfj,{kat) 



l/{l+2x)^-l 



This is faster than the Hubble rate for the scales fl3.2ip . as deduced earlier. 



(3.28) 



3.3 The long-string two-point function revisited 

We can now improve our earlier result for the two-point function found in 
Chapter 2 through the inclusion of gravitational radiation. Define 

/oo 
da e-'^'^iw+ia, t) ■ w+(0, t)) , (3.29) 
-00 
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where a is related to the actual length along the string by / = a(t)e(t)a. Recall 
that e{t) = a{t)~'^'" . Combining the stretching given by Eq. (]2.32p with the 
gravitational radiation found above gives 

^H{k, t) = --^aH{^, t) j^,H{k, t) . (3.30) 

at a at TQ^[K/ae) 

The derivative of the prefactor a{t)e(t) would give the first term with a plus sign 
but when the derivative acts on the correlator of the fluctuations the result is 
twice as large and negative. This equation is readily integrated, to give 

\nH{K,t) = -a\na{t) + \nf{K) - 6^(1 + + 2x)Gi2K^/^^''+^h'^^^+''^^^+^^^ . 

(3.31) 

The function /(/t) is obtained by matching onto the known result (13. ip at early 
times, giving /(/t) = ^c^/t~^~^^. Expressing the result in terms of the physical 
momentum scale k = K,/a{t)e{t) gives 

J{k,t) = H{K,t) = Ac^k~\kty^'' exp [-6^(1 + x)(l + 2x)G/i(A;t)i/(^+2x)] _ 

(3.32) 

Thus, unless we consider modes with such short wavelength that the expo- 
nential suppression factor kicks in, we have J{k) ~ /c~(^+^^). At length scales 
smaller than ~ (G/i)^+^^t this behavior is altered as a result of the inclusion of 
gravitational radiation, in such a way to make the two-point function smoother. 
In Fig. 13.11 we show the numerical Fourier transform of J{k,t), in terms of the 
correlation |37l EH] 

^ l-(p+(/,t)-p+(0,t)) _ (u;+(/,t)^)-(w+(/,t)-w+(0,t)) 
^-^OTT^ 2(1^) - 2(1^^ • ^^■^^> 



The smoothing at short distance seen in Fig. 13.11 resembles that found in the 
numerical simulations of Ref. [37], but here it is due to an additional physical 
effect, gravitational radiation, whereas in Ref. [3Z| such smoothing was found even 
without gravitational radiation]^ In Fig. 13.11 we have assumed that the enhanced 
straightness of the strings at the shortest scales seen in Figs. 12.21 and 12.31 is a 
transient effect, so that the two-point function follows a simple power law down 
to the gravitational radiation scale. 



^As an aside, the correlator in Fig. 13.11 goes as P at short distance, corresponding to an 
ensemble of smooth functions, whereas it appears to be close to in Ref. [37], corresponding to 
a random walk in p±. 
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Figure 3.1: The effect of the inclusion of gravitational radiation on the long-string 
two-point function is to smooth the string on scales below ~ 20{Gfiy~^'^^t. The 
solid line corresponds to the matter epoch and the dotted line to the radiation 
epoch. The value G/i = 10~^ has been adopted for illustrative purposes. 

The exponential falloff in the world-sheet two-point function at large wavenum- 
ber implies that the strings are very smooth at short distance. Inspection of 
Fig. 13.11 indicates that this smoothing sets in at 



in both the matter and radiation epochs. The exponent we obtain is smaller than 
the estimates in [68j: 1.2 in the radiation dominated era and 1.5 in the matter 
dominated era here, compared to 3/2 and 5/2 respectively. The difference is 
primarily from using a more accurate model of the effect of stretching, leading to 
a power spectrum with a slower falloff. The length /qw is the approximate size of 
the smallest loops. In particular, the divergent loop production found in Ref. [58] 
is cut off at this scale. 

The gravitational radiation back-reaction on the two-point function {p± ■ p±) 
obviously carries over to the two-point function for the fluctuations (w-|- ■ w-|-). 
These are fixed by the relations fl2.28p and (12.331) . or equivalently fl3.ip . but only 
up to an undetermined function 0('u) = |(w^(-u)): 



lew ~ 20(6-/1) 



(3.34) 



(w+(m) -w+Cu')) = -A 



u 




+ (j){u) + (f){u) , for \u — u \ ^ /gw • (3.35) 
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To specify the function one needs to choose how to spht p± into the slowly- 
varying P± and the fluctuation w±. A natural choice is to let the self-correlation 
of the fluctuations to be independent of the point on the string, i.e. is a constant. 
This constant is related to the behavior of the two-point function at large distances 
as can easily be seen by considering the Fourier transform of the above equation. 
The transform of the first part was essentially computed in (13. 8p . The constant 
piece (p can then be obtained by inserting a cutoff L representing the separation 
between the long-distance and short-distance parts of the configuration, taken 
around the correlation length: 

= ^r^ Re\{i\u-u'\+ Lf' -{i\u~u'\+lG^f''] , (3.36) 

and similarly for w_; recall that c^ = 2 sin(7rx)r(l + 2x). We have also introduced 
a smooth short- distance cutoff at the scale /gw ~ 2Q{GnY^'^H to take account 
of the effects of gravitational radiation, as found above. In the range /qw ^ 
\u — u'\ <^ L expression (I3.36P reduces to 

(w+(m) ■ w+{u')) ^ A (L'2>^ -\u- u'\^^)/t'^ , (3.37) 

with L'"^^ = L'^^/cosnx- So we see that our cutoff prescription corresponds to 
having (p = A{L' /tY^. This expression for the two-point function will be used in 
the next chapter to compute the rate of small loop production. 



3.4 Small scale structure and cusps 

Our work also allows us to address an old question, the interplay between the 
small scale structure and the cusps ^76j. A cusp is a point on a string that reaches 
momentarily the speed of light and where simultaneously the cx-parametrization 
becomes singular. We shall discuss this in more depth in Chapter 4 but we can 
anticipate that cusps arise when the functions p+(m) and P-(f) cross on the unit 
sphere [10]. The size of the cusp is the inverse of the 'velocity' of the functions p± 
when they cross, their rate of change with respect to the length along the string. 
Now, the average velocity of these functions on the long strings is (going again to 
the gauge e(t) = a{t) = 1) 

V = (9.p+(a,t) ■ 9.p+(a,t)) = J ^K'H{K,t) = d^{Gfi)-'^'-^^^'+'^h-' , 

(3.38) 
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where 

= Ac^^^mi - + 2x)] [6,(1 + + 2x)]^(^-^)(^+'^) (3.39) 

takes the value 0.006 in the matter era and 0.008 in the radiation era. 

On the long strings in the network, if the functions p+{u) and P-(f) were 
smooth and varied on the scale of the correlation length (which is not far be- 
low the horizon scale), then the typical velocity of any crossing would be the 
inverse correlation length, and this would be the size of cusps. What we see from 
Eq. (I3.39P is that such slow crossings and large cusps are impossible. Superim- 
posed on the slow motion is the irregular motion from the short distance structure, 
so all crossings occur at the much larger velocity V ~ 0.1{G jj,)^^^^^^^^^^^h^^ , and 
all cusps have a size of order 

~ 10(6-/2) , (3.40) 

which is very much smaller than the correlation length. Effectively each large 
cusp breaks up into a large number of small cusps: the high fractal dimension 
of p+(m) and P-(f) above the gravitational radiation length, noted in Ref. [58j, 
implies that these curves will cross many times near the would-be large cusp. We 
can estimate the number: the total length of each curve p+(m), P-(f ) during one 
Hubble time is Vt ~ 0.1(G'/i)-(i-'^)(^+2x\ The total number of crossings is roughly 
the total length of either path divided by the solid angle An of the sphere, and so 
up to numerical factors it is large to the same extent that the size of each cusp is 
small. 
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Chapter 4 
Loop formation 



The model presented in Chapter 2 assumes that stretching is the dominant 
mechanism governing the evolution of the small scale structure. Recall that the 
output of this model is a power law describing the approach of the two-point 
function to straightness as we consider smaller separations along the string. How- 
ever, the simulations of Ref. [37J show strings even straighter than our stretching 
model suggests, at the smaller scales. We have speculated in Section 12.41 that this 
discrepancy could be due to loop production. Given the smoothness of the strings 
on short distances, one might have thought that loop production on small scales 
would be suppressed, but we will see that this is not the case. 

In this chapter we include the production of small loops in our model. We 
will consider the stretching model as a leading approximation and add in the loop 
production as a perturbation. Loop reconnection will still be ignored, based on 
the standard argument that this is rare for small loops [8|. 

We shall begin by reviewing some useful results regarding loops, cusps and 
kinks in Section 14. 1[ We then proceed to compute the rate of loop production. 
This calculation is contained in Section lT2l Section H73l is devoted to the discussion 
of the results and includes some speculative comments as well. 

The treatment presented here is based on Refs. ^Hll62], with Joseph Polchinski 
and Florian Dubath. 

4.1 General results: loops, cusps and kinks 

Loops much smaller than the horizon size evolve in an essentially fiat spacetime 
so we are able to use null coordinates u = t + a and v = t — a. In this case the 
left- and right-moving vectors p± depend only on m or f , respectively. They are 
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related to the spatial embedding of the string, x(t, a), simply by 

P+(m) = 29„x(M,f) , P-(f) = 2d^:s.{u,v) . (4.1) 

Recall that these vectors are constrained to live on the unit sphere. 

It is easy to derive the condition under which a loop will form by a self- 
intersecting long string. Assuming that the reconnection probability is P = 1, 
a loop of length I will form if there exists a worldsheet point (t, a) for whichu 
x(t, a — 1/2) = x(t, a + 1/2), or equivalently, 

r<T+l/2 

/ da'd^'x{t,a') = . (4.2) 

Ja-l/2 

Thus, using = du — and the relations (14. ip we conclude that the condition 
for the formation of a loop of length / is 

L+(m,/) = L_iv,l) , (4.3) 

L+(?i, /) = / du' 'p^{u') , L_(t>,/)= / dv' 'p-{v') . 

Ju-l/2 Jv-l/2 

Stated differently, the unit vectors p-t must average to the same value in order to 
yield a self-intersection. 

Note also that a loop of length I actually has period T = 1/2 in time. The 
simplest way to see this is to recall that in fiat space signals propagate in both 
directions at the speed of light. Left- and right-moving signals originated at some 
point (Jo on the loop will meet again at the 'opposite end', ctq + 1/2. The period 
is equal to the distance traveled, which yields half the length of the loop. 

Let us now consider the particular case when the two curves p± cross. Using 
(9tX = |[p+ + p_] it immediately follows that such points on a string (called cusps) 
move (instantaneously) with the speed of light [79J, (9jx)^ = L Correspondingly, 
the tangent along the string 9o-x momentarily vanishes at the tip of the cusp, 
signaling the singular behavior of the parametrization of the string by cr. The 
highly relativistic velocities around these points translate into a strong emission of 
gravitational waves along the direction of motion of the cusp [271 180] . This feature 
has been used by those authors to study the detectability of such gravitational 
waves. In Chapter 6 we will consider instead gravitational waves originating from 
the lowest harmonics of small (and relativistic) loops. 



^Here we choose to center the segment of the long string that originates the loop around 
(T = 0. 
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When two smooth segments of string intersect and reconnect this resuhs in 
two kinked segments, as shown in Fig. 12.11 A kink is a discontinuity of the 
tangent vector c^o-x. In terms of the left- and right-moving unit vectors this is 
identified by a discontinuous jump in both p±. As long as the approximation of a 
flat spacetime can be maintained, each discontinuity evolves into a pair of kinks 
traveling at the speed of light in opposite directions along the string. A simple 
way to understand why this happens is to note that one of the kinks is associated 
with a discontinuity in p+(m) and so is left-moving whereas its 'twin' kink is 
related to a discontinuity in P-{v) and therefore is right-moving. More physically, 
such evolution shortens the length of a string between two points surrounding the 
initial kink and consequently reduces its tension. 

4.2 The rate of loop production 

In this section we aim to determine the rate of loop production, which occurs 
when a string self-intersects. More specifically, we want to compute the number 
of loops df/ with invariant length between / and / + dl originating from a self- 
intersection at a string coordinate between a and a + da during a time interval 
{t,t + dt). 

4.2.1 Preliminary considerations 

The calculation of the loop production rate requires separating the p± into 
a long-distance piece and a short-distance piece. In Ref. [5H] the long-distance 
part was simply averaged over the unit sphere, with a weight factor chosen to 
give the correct mean value of p+ ■ p_. Assuming the stretching model derived 
in Chapter 2 remains valid down to arbitrarily small scales, it was found that the 
loop production function diverged at short distance. Clearly, this divergence must 
be cut off because the total rate at which string length goes into loops is fixed by 
energy conservation. In the approach of Ref. [58], this condition saturated at a 
length scale just one order of magnitude below the horizon. However, comparing 
with the numerical simulations of Refs. [l5l [38] that treatment did seem to find 
the correct functional form for the loop production function. 

An essential improvement was made in Ref. [62], leading to a very different, 
and simpler, picture. In that work the long-distance part was given a classical 
(m, f )-dependence, and then averaged over the short distance part with a weighting 
corresponding to Eq. (12.361) . The inclusion of the time-dependence reveals that 
loop formation leads to a strong correction to the distribution of long-distance 
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Figure 4.1: The functions p± are written as a fixed long-distance piece (bold) 
plus a random short- distance part. The classical cusp has been defined to lie at 
{u,v) = {0,0). 

configurations, which was not accounted for in the earlier approach |58jl^ 

A typical configuration for the p± vectors is shown in Fig. 14. 1[ The long- 
distance configuration is in bold, and the total configuration is dashed. Recall 
from Section [27il that . while the string configuration approaches fractal dimension 
1 at short distance, the graph of its tangent has large fractal dimension 1/x- In the 
figure, we have shown a long distance configuration in which the curves p+{u) and 
p_.(f) cross (at M = and v = 0), meaning that there is a cusp in the spacetime 
evolution of the string. It is evident that the effect of the short distance structure 
is to turn this one large cusp into many small cusps. Therefore we see that the 
formation of small loops is closely associated with the cusps. 

When condition (14. 3 p is satisfied the segment of length /, centered at {u,v), 
breaks off as a loop. For small / the functions p-t will have small variation along 
the segment, so small loops will form where these are approximately equal [HUES]- 
Following this discussion, it is clear that at early times, when u and/or v are 
sufficiently negative, the curves will be outside the cusp region and the production 
of small loops will be negligible. As the cusp begins to form, production of small 
loops begins to increase. We will find that loops of all scales begin to form at the 
same time, extending all the way down to the gravitational radiation cutoff. Thus 
we are led to conclude that the production of loops down to the gravitational 

^ These same issues were faced in the pioneering work |52) . but our framework differs. 
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radiation scale is inevitable, and not a transient effect. 

The loop production cuts off when the probability for a given bit of string to 
be incorporated in a loop reaches 1. We will see that this occurs long before the 
classical cusp point (u, v) = (0, 0) is reached. In the approach of [58j the region 
where the long-distance parts of p+(u) and P-(f ) were parallel gave the dominant 
contribution, but in a correct treatment it is simply absent. 

In the following section we will carry out the calculation just described. As 
further support for our picture, we will obtain a prediction for the small loop 
production function that fits well with the simulations by two groups. The fit 
actually works better than anticipated, as there is no sign that it is distorted by 
further fragmentation; we speculate as to why this might be. 

4.2.2 Calculation of the rate 

To perform the calculation of the rate of loop production from long strings, 
we separate the left- and right-moving unit tangent vectors into their long- and 
short- distance parts, 

P± = 1± + s± . (4.4) 

Since the production of small loops takes place near the cusp, we take the simplest 
cusp form for the long-distance structure, 

\+{u) = z{l-V^u^y^'^ + Y+u^z + V+u-^V^u\ 

= z(l- VV)^/^ + V_f ^z + V_t;- , (4.5) 

with V_(_ ■ z = V_ ■ z = 0. We have chosen a coordinate frame in which the motion 
of the cusp is directed along the z-axis and we are expanding to order u^, . 
The slopes V± have units of inverse length, and are reciprocal to the size of the 
cusp. They will be of order the correlation length, somewhat less than the horizon 
length. 

The vectors 1± encoding the long-distance structure play the role of the vectors 
P± defined in the small fluctuation approximation of Section 12.3. 1[ In particu- 
lar, they are normalized to unity. The short-distance structure parts are thus 
expressed in terms of the small fluctuations as 

s± = w± - . (4.6) 
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Since the fluctuations must be orthogonal to z, we can insert expressions (14.51) 
and ( 14. 6p into Eq. ( 14. 4p and obtain, again to quadratic order, 

1. 2 

p± = z + y± - ^zy± , 

y+(n) = V+n + w+(n), y_(i;) = V_i; + w_(i;) . (4.7) 

For the fluctuation two-point function we shall use (I3.37p . Actually, the classical 
(m, f )-dependence given to the long-distance parts of p± should add a term of the 
form \V^{u — u'Y to the right-hand side. However, this is quadratic and we can 
ignore it when considering the production of loops as a result of the small scale 
structure. 

The number of loops with length between / and / + dl formed by self-inter- 
sections of a long string occurring within an infinitesimal worldsheet area dudv is 
given by 

dJ\f= {5^^\'L+{u,l) -L4v,l)) \det J \) dudv dl , (4.8) 
where J is the Jacobian for the transformation {u, v, I) L_|_ — L_: 

p+iu + l/2)-p+iu-l/2) 
J{u,vJ)= p-{v-l/2)-p^{v + l/2) 

|[p+(n + 1/2) + p+iu - 1/2) -p.{v + 1/2) -p.{v- 1/2)] 

(4.9) 

This formalism is as in Refs. [531 152] . 

In |5H| the computation proceeded with an estimation of the expectation value 
of the product (14. 8 p by the product of expectation values, but here there are strong 
correlations: when the delta-function is nonzero, the third row of the Jacobian is 
much smaller than its mean value. We therefore go to new variables, separating 
w-i- into a piece constant on the segment and a piece with zero average on the 
segment, 

w+(?i') = W+ + a;+(n'), W+ = - / du'w+{u'), 

1 rv+l/2 

= W_+a;_(t;'), W_ = - / dv'w^v'). (4.10) 

' Jv-l/2 



The variables u}± and W-t depend on the parameters u, v, I of the loop, but we 
leave this implicit. We will be interested in the loop production in the range 
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^GW I L, SO we use the form fl3.37p . One then finds the two-point functions 



t2x^W+ ■ W+) 



-A\u' - u"\^^ + f{u') + f{u") 



where 



A 



(2x + i)/ L 



{1/2 + u' - + {1/2 -u' + uf^^^ 



/2X+1 

2x + 2, 



(4.11) 
(4.12) 
(4.13) 



(4.14) 



is defined only for points on the segment forming the loop, i.e. u — ^<u'<u + ^. 

Now let us express the rate in terms of these quantities. First, for the transverse 
parts of L-t we have 

= /(W+ + V+m) , = /(W_ + V_f) ; (4.15) 

the transverse part of the delta-function sets these equal L:[: = L:*: = L-*-, and we 
use this in evaluating the remaining terms. For the z component. 



2 
1 



/ - — L 

21 



u+l/2 

u-l/2 
2 



1 ^+1/2 
^ Ju-l/2 



U 



(4.16) 



where we have used the fact that the mean value of a;+ is zero. In the integrand, 
the fiuctuations of a;+ are of order l^, while the classical term is of order /, so we 
can drop the latter for small loops. Thus, 



u+l/2 



du' [u;+{u')f 



1 



u-l/2 



v+l/2 



dv'[iv4v')f = 0(/i+2^) . (4.17) 



v-l/2 



This takes care of the delta function in (14. 8p . Now we turn our attention to 
the Jacobian. For the transverse part of the first row of J, 



pi(M + //2)-pi(M-//2) 



uj+{u + 1/2)- u)+{u -1/2) + V+l 

u;+(m + 1/2) - uj+{u - 1/2) + 0{l) . (4.18) 



Similarly in the second and third rows, in the transverse terms we can replace p± 
with uj±, after imposing L:[: = L:*:. In the z component of the first row, 

pI{u + 1/2)-pI{u-1/2) = -jL^ ■[uj+{u + l/2)-Lj+{u-l/2)] 

-\ (K(t. + l/2)f - K(^ - l/2)f) + 0{l) . (4.19) 
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The first term is of order and tlie second of order l"^^, but the first actually drops 
out. In all three rows of J one finds the same pattern, = L"*" ■ + Thus 
the first term in the z-column is linearly dependent on the other two columns, 
and drops out of the determinant, leaving with replaced by 

We have now expressed the terms multiplying 5^'^\Ij-^ — L:*:) all in terms of lj. 
The cross-correlation fl4.12p is of order l"^^, smaller than the geometric mean L'^l^ 
of the diagonal correlators fl4.11ll4rT3|) . We can therefore ignore it, giving 

(5(=^)(L+ -L_) |det J|) = (5(')(L^ -L^)) X -LI) |det J|)p^_^ . (4.20) 

Thus we have achieved our aim of factorizing the expectation value. In the sec- 
ond term the delta-function scales as inversely to its argument, while the 
columns of J scale as l^, l^, /^^, giving the overall scaling, as found in |58j . 
This second term is evaluated in Appendix A with the assumption that the short- 
distance structure is Gaussian. One finds that 

{6{Ll - U) |det J|)px^^ = r^- j , (4.21) 

where rj 2.3 in the radiation era and rj ^ 1.3 in the matter era. As we will 
discuss in Chapter 6, the Gaussian approximation may not be valid here, but this 
does not alter the scaling with /. 



It remains to evaluate 



<5^(Li-L^)> = j 



giq.(Li-Li) 



(27r)2 
(2vr)2 



j-2Xp-(V+«-V_i;)2t2x/2^L'2x 

. (4.22) 



In all. 



dN = CL'-2Xe-(v+"-v-)^*^V2^L'2x dudvj^, (4.23) 

where (with the Gaussian approximation) C = 0.35 in the radiation era and 
C = 0.2 in the matter era. These constants are actually independent of A so 
the normalization of the two-point function only enters the calculation of loop 
production through the exponential suppression factor. 



^Note this implies that is replaced by 1 — ^(j^±, since the vectors p± are normalized to 
unity. Furthermore, the constant term may be dropped because it cancels out in all calculations. 
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The result above indicates that the total rate of string loss J I dM diverges 
at small / for x ^ 0-5, as is the case in both the radiation and matter eras. 
Thus, loop formation can have a significant effect on the evolution of small scale 
structure and could invalidate the assumptions of our model but we will argue 
against this in Section H. 3. 2[ Of course, the total rate of long string lost into loops 
is bounded. Indeed, the small-/ divergence must be cut off at the gravitational 
radiation smoothing scale found in Chapter 3. But, for now, let us verify the 
statement made in Section l4.2.1l that segments on the long string with the potential 
to form cusps are carried away by loops before the cusp actually forms. 



4.2.3 The normalization 

Consider a left-moving point with given u. The total probability per unit v 
that this point be incorporated into a loop is 



^ = CL'-2^e-(^+"-^-^)'*''^/2-^^"^ j ^ . (4.24) 



Note that we have replaced J du I to count the loops containing the given 
point. This diverges at small / because we have not yet taken into account the 
smoothing due to gravitational radiation in the two-point function fl3.36p . The 
smoothed two-point function is quadratic at u' — u ^ 0, corresponding to the 
form fl3.37p at x = 1: at the shortest distances the divergence is gone. Thus we 
cut off the integral at /gw to get 

dv 

At early times, where v is large, the probability of a loop containing u is small. 
However, Iq^'^^ is large compared to the other dimensionful quantities, so as we 
integrate in v we soon reach total probability 1 for a range of values of u near 
the cusp: this portion of the string is removed by loop production (see Fig. 14. 2p . 
Therefore, for each value of the w-coordinate, the maximum value of v that may 
still not be included in a loop is the solution (when it exists) of the equation 
V{vm.ax) = 1, where 

r(viu)) = ( ^ ) / dv e-(v+-V-.)^*^x/2^,..2. _ ^^ 26) 



:i - 2x)Igw V L' 

We see that the excised region depends on the the angle 6 between V_|_ and V_ 
for each cusp, not to mention the ubiquitous G^. To obtain the normalization of 
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Figure 4.2: The shaded region in the (-u, f)-plane is excised by loop formation. 
The cusp occurs aX u = v = and the curve that dehneates the excised region is 
Vmax{u)- Both u and v are expressed in units of the cosmological time t. For this 
particular example the values G/i = 10~^ and cos(^) = 2/tt were used, where 9 
denotes the angle between V+ and V_. 

the loop production function would require the exact knowledge of the probability 
distribution for the angles 6 as well as the density of cusps. One can certainly 
obtain normalizations in accordance with the simulations but there are too many 
free parameters. 

Thus, it is beyond our present capabilities to determine the normalization 
of the loop production function directly from the knowledge of the two-point 
function. Nevertheless, one can still make some progress by employing energy 
conservation. The total rate of string loss is bounded: the long string network 
cannot loose more energy than it initially contains! In a scaling solution, the total 
amount of long string in a comoving volume scales as ioo oc a^/t"^, while stretching 
alone would give ioo (x a". The rate of string lost to loops must be 



Averaging expression fl4.23p over the ensemble of segments and noting that dudv = 




(4.27) 
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2dtda we must have an average rate of loop production of the form 

W c //^'^-' 



dtdadl \t 



(4.28) 



We were not able to compute the normalization c from our stretching model but 
now we can fix it by equating the rate of string length going into loops to the rate 
of length lost by the long strings. 



(4.29) 



9t /loops J low dtdadl Ji^^ t3 yt 

C / Lew 



We have inserted a cutoff at the gravitational radiation scale low- Equating the 
two rates gives 

c=(l-2x)[2-2Kl+t)^)] (^^j . (4.30) 

So we conclude that it is the gravitational radiation (or at least some UV cutoff 
on the sizes of loops produced) that determines the normalization of the loop 
production function. 



4.3 Discussion 

These are the main results of this chapter. First, loops of all sizes down to 
the UV cutoff form simultaneously, rather than in a cascade of fragmentation. 
Therefore the production of the small loops is a robust physical result. The 
distribution of the loop lengths obeys a simple power law with the exponent 
determined by x- 

The production of very small loops takes place when p_|_ ~ p_: that is, near 
cusps on the long strings as proposed in Refs. [81], [82]. Indeed, all the loops that 
we consider are rather smaller than the correlation length, so the functions p+ 
and p_ are each somewhat localized on the unit sphere, and necessarily in the 
same region since L+ = L_. In this sense all of our loops are produced near cusps. 

The integral of the probability (I4.25P for a point on the string to break off 
becomes large long before the cusp u = is reached. In Ref. [58] this fact was 
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not accounted for and this resulted in a large overcounting of loops formed by 
self- intersections. 

It is possible that there is a second population of loops that form at much 
longer scales and we will speculate on this issue below. But first let us compare 
our results with some of the most recent simulations. 



4.3.1 Comparison with simulations 

The results that we have found have some notable agreements, and disagree- 
ments, with simulations. One success is an apparent agreement with the simula- 
tions of Ref. 08] for the distribution of loop sizes: that reference finds a number 
density of loops per volume and per length dn/dl oc with p = 2.41 in the mat- 
ter era and p = 2.60 in the radiation era. We have exponents 3 — 2% = 2.5 and 2.8 
respectively. Our exponent is for the production rate rather than the density, but 
in this regime the density is dominated by recently produced loops and so these 
are the same. 

Let us verify this, and also obtain the relative normalizations for the two 
quantities. The number of loops per comoving volume of length between / and 
I + dl is unaffected by the expansion of the universe, and we are assuming that 
we are at scales where gravitational radiation can be neglected. The number then 
changes only due to production: 

d { .dn\ a!^ c fl\^^ /. o.\ 

37 « 37 = 7^ « 7 • (4.31) 



dt\ dl J 72^2 p y^ 

(Recall that 7 is a dimensionless constant: the scaling value of the long string 
length per unit volume is (7t)~2). This integral is dominated by recent times as 
long as 3z/ — 1 — 2% > 0, as holds in both eras. Integrating over time then gives 

dn c fl\'^^'^^ 

an c IL\ ^^^^2) 



dl (3z/ - 1 - 2^)72^^ \t 
In the notation of Ref. [3^ the normalization is defined as 

Co = 7 r3(l - u)-^-^'' , (4.33) 

(3z/-l-2x)7'^ ' ^ ' 

where the last factor comes from the conversion from dn to t. 

In Section l4.2.3l we determined the normalization c in terms of the gravitational 
radiation cutoff. Ref. [38] does not explicitly consider gravitational radiation but 
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their simulations have a minimum resolution: loops below a certain length are 
artificially removed from the network. As we will see in the next chapter, this 
effectively introduces a gravitational radiation scale consistent with the values 
low — 1.2 X 10~H and lew — 6 x 10~^t for the matter and radiation eras, re- 
spectively. Using 7 = 0.59, 0.30 respectively in the matter and radiation eras p?7] . 
this gives Co = 0.07, 0.03. These are smaller than found in Ref. [38] by factors of 
1.3 and 7. Given the crudeness of our estimate, the disagreement by less than an 
order of magnitude is a pleasant surprise. 

A different group [46j has performed simulations that, due to a doubling trick, 
have been able to reach a larger dynamical range than in previous numerical works. 
Their results relate more directly to ours since they consider explicitly the rate 
of loop production. Those authors express this in terms of where x = l/t 

and the loop production function (LPF) f{x) is related to the rate through 

m = ■ (4-34) 

These simulations show a power law distribution above the UV cutoff, with ex- 
ponents that match rather well with our model in both the matter and radiation 
eras. The issue of the normalization is a more speculative story. First, the nor- 
malization of the loop production function in |36] is still decreasing when the 
simulations stop running. Second, these simulations do not include gravitational 
radiation and have no minimum resolution, so according to Eq. (14.301) it seems 
that the normalization should approach zero at late times. Indeed, again by con- 
servation of energy, if the LPF retains the power law behavior down to smaller 
and smaller loop sizes as the simulation progresses then its normalization must 
decrease. Since the smoothing by gravitational radiation is a physical effect, the 
evolution of the LPF must stop when the scale Iqw is reached and in that case the 
normalization should be fixed by Eq. (14.301) . On top of the power law distribution 
we have been discussing, Ref. |16] also detects a peak at quite large loop lengths, 
I ~ O.lt and we shall comment on this below. 

As we have already mentioned, Ref. [33] also finds power law distributions 
for x^f{x). In Table 14.11 we compare the values of the exponents obtained in 
the simulations by both groups with the exponents predicted by our model. The 
simulations [38] and [16] thus agree reasonably well with each other and our model 
lies in between. 
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Ref. 06] 


Ref. [38] 


our model 


Radiation dominated 


^-0.95 


^-0.60 


^-0.8 


Matter dominated 


^-0.55 


^-0.41 


^,-0.5 



Table 4.1: Comparison of the exponents in the loop production function power law 
obtained from simulations by two groups with those predicted by our stretching 
model. 

4.3.2 Effect of loop production on the two-point function 

The model developed in Chapter 2 ignored the production of small loops as 
being a sub-dominant effect to stretching. However, we have found that, although 
the fractal dimension of the long strings approaches 1 at short scales, the rate 
of loop production becomes relevant in cusp regions. Does this invalidate the 
assumptions of our model? 

Recall there was a discrepancy between the two-point function as dictated by 
our model and the results of [37]. At the shortest scales in Figs. 12.21 and 12.31 the 
slope of the curve approaches and possibly even exceeds unity, corresponding to 
the critical value x = 0.5: if the two-point function is of this form then the small- 
loop production converges. Thus it is appealing to assume that it is feedback from 
the production of short loops that accounts for the break in the curved 

We believe this is not the case. If the loop formation were distributed on the 
long strings, there would be at least one additional effect, where the string shortens 
due to loop emission and more distant (and therefore less correlated) segments 
are brought closer together. However, the loop production is not distributed: it 
occurs when the functions p±(cr), as they wander on the unit sphere, come close 
together. Our picture is that whole segments in cusp regions are then excised. 
Thus, most of the segments that remain on a long string at a given time would 
have been little affected by loop production. 

Consider a very short segment on a long string. With time, its size as a fraction 
of the horizon size decreases, so it moves to the left on Figs. 12.21 and 12.31 If it 
experiences only stretching, it will follow the slope of the dashed line. Thus we 
would expect the two-point function to be approximately given by the stretching 
power law down to arbitrarily small scales, until gravitational radiation entersO 

^Note also that x = 0-5 corresponds to the functions p± being random walks on the unit 
sphere. This suggests this value might be some dynamical fixed point. 

^One effect that may increase x slightly is that segments that do not fall into the loop 
production 'hole' at p+ p_ may experience on average a larger value of a = — p+ • P- and 
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The loop production must have a strong bias toward removing segments with large 
fluctuations to account for the simulations. 

There is a somewhat related puzzle regarding the loop velocities, first noted 
in [58]. For a loop of length /, the mean velocity is 

V = ^(L+(n, /) + L_{v, I)) = jL+{u, I) . (4.35) 



Then 



9 A 

[l/tf^ . (4.36) 



(2x+l)(2x + 2 

Looking for example at loops with I = lO^^t, we obtain a typical velocity 0.985 in 
the matter era and 0.90 in the radiation era. These are significantly larger than are 
generally expected; Ref. [S3] gives values around 0.75 and 0.81 for loops of this size. 
We emphasize that this is not a consequence of our dynamical assumptions, but 
can be seen directly by using the two-point functions of Ref. [37J (the discrepancy 
is then even greater in the radiation era). Assuming that both simulations are 
correct, we must conclude that the two-point function on loops is very different 
from that on long strings, in fact less correlated. One possible explanation is that 
loop formation is biased in this way, but this could also be an indication of the 
complicated nature of the fragmentation: the final non-self-intersecting loops must 
contain segments which began on the long string many times further separated 
than the final loop size. 



4.3.3 Fragmentation 

The formation of loops from self-intersections on long strings is just the begin- 
ning of the story. After the primary loops are chopped off (and even before this 
happens) the process of loop production continues on smaller scales, unaffected 
at least initially by the separation of the loop from the infinite string. Thus, the 
primary loops fragment into smaller daughter loops. 

We have found that fragmentation is not necessary to produce the small loops, 
but it still raises a puzzle. Given the small-/ divergence of the rate of loop for- 
mation in cusp regions one might have expected a rather dramatic fragmentation 

so have a reduced two-point function — see Eq. (I2.20|) . However, it seems impossible that this 
could account for the factor of 5 discrepancy of slopes evident in Fig. 12.21 
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process. Indeed, if we were to apply our same calculation to study fragmentation 
of these small loops, we would find again a strong tendency to fragment into loops 
at the gravitational radiation scale, giving a sharp peak there rather than a power 
law distribution. 

In fact, we believe this is not precisely the situation. Rather, a non-self- 
intersecting loop will occasionally form and this length of loop will then be lost 
to the fragmentation process on smaller scales. The point is that the short dis- 
tance structure on these small loops is not that typical of the long strings, but is 
dominated by two or more large kinks moving in each direction. One kink in each 
direction forms with the loop, and the remainder are present in the pre-existing 
distribution. From Fig. 14.11 we expect that if there are any cusps on these loops 
then the regions containing them will dissolve into small loops. However, if there 
are parts of the loop where p+ is sufficiently far from the p_ curve, and vice versa, 
then these will survive. Thus we expect that if the functions p± are plotted, there 
will be significant gaps in each, corresponding to kinks which form simultaneously 
with intercommuting events, and the surviving curves will avoid each other. 

It was noted in Ref. |58] that the tails of the small scale structure distribution 
are not Gaussian but are dominated by a single large kink (again, this does not 
alter the scaling with /). Thus these kinks are likely associated with the tail 
of the loop production function. If so, the functions p± on the small loops are 
dominated by large jumps, from the kinks, and will be unlikely to have cusps 
where smaller loops can form. This conforms to the conclusions of [H], where 
numerical simulations of the fragmentation of loops was studied. There it was 
found that there is an end to the fragmentation process, meaning that all loops 
eventually reach non-intersecting trajectoriesJ§ 

Finally, the fragmentation need not change the scaling of the LPF as long as 
the long-string correlation function remains a power law, because the shapes of 
the produced loops, and the resulting fragmentation, all scale. Thus, we do not 
consider the agreement between the exponents characterizing the loop production 
obtained from our model and from the simulations to be accidental. 



4.3.4 Large loops 

We have already mentioned that Refs. HE] find two populations of loops, 
one of which scales at a = l/t ~ 0.1, and a second at much shorter scales which 
is still evolving toward smaller lengths. We believe the latter corresponds to 
the power law distribution borne out of our small scale structure model. The 

^However, we note that the spectrum of perturbations used to define the initial loops in this 
study does not resemble our distribution for cosmic strings in expanding spacetimes. 
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production of large loops is outside our model, but these also seem inevitable. In 
the course of the evolution of the network it must happen at least occasionally 
that a large loop will form and survive reconnection. We must leave it to the 
simulations to determine the relative weight and average size of these loops since at 
these large scales the cosmic string evolution is highly non-linear, but the argument 
given in the previous section explains why they do not fragment entirely into much 
smaller loops. An account of these large loops has been given in Ref. [85] . 

It seems that the interpretation of the small loops as being transients is im- 
plausible. First, there is no sign in Ref. that the total fraction of string length 
going into the small loops is diminishing in time. That was the case for Ref. [15] 
but those simulations were done in flat spacetime. Second, if this were a tran- 
sient effect then its natural time scale should be small in proportion to the size 
of the loops. Instead, these loops are still being produced after several Hubble 
times, showing that some longer-distance process is continuing to feed the small 
loop production. Indeed, this mechanism has already been identified: the small 
scale structure on long strings and large loops, which originates at the horizon 
length and then is carried down to smaller relative lengths by the expansion of 
the universe, leads to production of loops on arbitrarily small scales. 

Thus we are led to suggest the following picture. After a large loop is formed, 
production of smaller loops continues near its cusp regions. In the end, there is a 
large loop without cusps or self-intersections, reduced in size by the excision of the 
cusp regions. There is also a population of much smaller loops; based on the non- 
scaling seen thus far [45l|46], we conjecture that these will be at the gravitational 
radiation lengthQ Inspection of Fig. 4 of Ref. [16] suggests that as much as 80% 
of the total string length goes into the small loops, with 20% remaining in the 
large loop. 

In conclusion, it appears that numerical and analytic methods are converging 
on a firm picture of the loop production, with one peak near the horizon scale 
and one near the gravitational radiation scale, with perhaps the larger fraction of 
string in small loopsj^ 



^In the absence of gravitational wave smoothing, these could even be at the scale set by the 
thickness of the string [74] , 

®We are following Refs. [TSKSS] in putting the UV cutoff at the gravitational radiation scale. 
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Chapter 5 

Properties of the small loop 
population 

As mentioned in the Introduction, the large scale properties of cosmic string 
networks are well understood from simulations [TUl El [H] and there are analytic 
models that accurately capture those features [551 ES]. However, this is not the 
case regarding the issue of small scale structure. For example, the existence or 
not of a scaling regime for the smaller scales (for example, the typical size of the 
small loops) still remains an open question [HSJ. 

In this chapter we collect relevant properties of small loops and attempt to 
characterize such a population. The phenomenology of small loops is addressed 
in Section 15. 1[ In Section 15.21 we consider the number density of loops - it is the 
small scale structure that is responsible for their production and so the number 
density provides a measure of the former - and demonstrate that this quantity 
scales, both neglecting and including the effects of gravitational radiation. 

Most of the material contained in this chapter is drawn from Refs. [59l l60l [6T]. 
with Joseph Polchinski and Florian Dubath. 

5.1 Phenomenology of small loops 

5.1.1 Gravitational radiation from small loops and their 
lifetimes 

Let us start by briefly reviewing the issue of gravitational radiation emitted by 
small cosmic string loops. For the "vanilla" cosmic strings, the rate at which the 
loops radiate energy in gravitational waves determines their lifetime and this will 
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be used in Chapter 6. In addition, when the radiation from the whole population 
of loops is added up, it gives rise to a stochastic gravitational radiation back- 
ground [86] and this can lead to detection or better bounds on the dimensionless 
string tension, namely by accurate measurements of pulsar timing [281 EZl [33]. 

A heuristic derivation of the gravitational radiation power can be obtained by 
using the quadrupole formula [86] . 

^^G(^^y ^GmW, (5.1) 



dt \dt^ J 

where / ~ mi'^ is the quadrupole moment, m ~ yu£ is the mass of the loop and 
uj ~ is the characteristic frequency for a loop of size i. Therefore, the power 
radiated in gravitational waves is expected to be roughly 

P = TGfi^ . (5.2) 

Thus, the rate of energy radiated by small loops is proportional to Gfi"^, just like 
for the long strings. The proportionality coefficient F depends on the trajectory 
of the loop in question but not on its overall size. 

The validity of this heuristic argument can be questioned since the quadrupole 
formula ceases to be accurate in the presence of relativistic motion of the source, 
which is the typical case for oscillating small loops. However, more careful treat- 
ments confirm the result ( 15. 2p . The values suggested by numerical [HHl [89] and 
analytic [90l |9l] studies are T ~ 50. 

Since a loop of length / has energy fil, equation (15.21) implies that loops shrink 
at a rate 

ft = ' ^^-^^ 

and so a loop which formed with an initial length k has a lifetime given by r = 
{rGfi)~^li, if gravitational radiation is the principal decay mode. It is customary 
to parameterize the size of small loops as I = eTGfit. Thus, the parameter e 
corresponds to the lifetime of the loop in units of the FRW time. Given the power 
law distribution for the small loop production found in Chapter 4 we expect a large 
population of loops accumulating at the gravitational radiation scale. Therefore, 
recalling equation (13.341) from Chapter 3, we have found 

e ^ 0.4(G/i)°'^ in the matter era , (5.4) 
e ~ 0.4(Gyu)*''^ in the radiation era . (5.5) 

For loops containing ideal cusps the angular distribution of the power radiated 
is divergent along the direction of motion of the cusp [8^. Nevertheless, the 
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divergence is integrable and the total power remains safely finite. Gravitational 
radiation is expected to round off the tip of the cusp on very short scales and 
so it never really reaches the speed of light but comes quite close. The angular 
distribution of the power radiated is then bounded from above but this effect can 
still lead to very strong bursts of gravitational waves, which can be observed, 
in principle, at GW detectors [801 ESI [93]. A similar idea will be exploited in 
Chapter 6, but instead of the high frequencies typical of cusps we shall consider 
the lowest harmonics of the loops which get Lorentz boosted by the high velocities 
of small loops, a subject to which we now turn our attention. 

5.1.2 Loop velocities and cusps 

For the phenomenology of the small loops there is another important property, 
first noted in Ref. |58]: they move extremely fast. The point is that the functions 
p+{u) and P-(f) must have equal mean values on the loop (the condition for the 
loop to form). If the loop is small, so that these functions have limited range, 
this implies that both remain near one point on the sphere, and the velocity (the 
average of p+{u) and P-(f) over the loop) will be very close to unity. This is 
easy to understand since the whole loop is formed from a small region around a 
cusp and such points approach the speed of light. From Eq. (14.361) . the Lorentz 
contraction factor is 

7 = (/A)-^[(l + x)(l + 2x)M]^/^ (5.6) 

This is ~ 0.8(G'/i)~'^'^'^ in the matter era, and is of order 10^ for tensions of interest; 
in the radiation era it is ~ 1.0(G/i)~°'^^ and is of order 10^. Since this factor is so 
large, we must be careful to recall that the specific definition of length is 

dl = a{t)e{t)da = a{t) - ±'^)da (5.7) 

and so / is 7 times the rest frame length l^^st- It is equal to twice the inverse period 
of the loop in the FRW rest frame. 

To see one effect of the boost, in Ref. [94j the total stochastic gravitational 
wave spectrum from small loops was considered. The frequency spectrum simply 
scales inversely as the size of the loops if they are at rest, uj oc l//rest- In the FRW 
rest coordinate this frequency is reduced by a factor of 7. However, the large 
boost pushes essentially all of the radiation into the forward direction, where the 
frequency is enhanced by a factor {l—v)~^ ~ 27^. Thus the gravitational radiation 
spectrum corresponds to an effective loop size 

SeS ~ 7"'^rcst/2 ^ 7~'^/2 • (5.8) 
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This is ^eflf ~ 0.3(G/i)^'^^ in the matter era, meaning that Zcfr ~ 15(G'/i)^'^^t. In 
the radiation era Ees ~ 0.2(G/x)°-'^'', meaning that Ics ~ lO{Gfiy'^H. We should 
note that Ref. [M] considered the sensitivity in the whole [Gfi, e) plane, and we 
are simply highlighting one line in this plane. It is notable that Advanced LIGO 
is more sensitive than LISA to these smallest loops, because their small size puts 
them above the LISA frequency range. In fact, for interesting values of the tension 
the loops will produce a periodic signal in the LIGO frequency band, which may 
be observable at Advanced LIGO. Such a study has been conducted in [BDJ and 
we reproduce it in Section 16. 2[ 

Recall further that the cusp region is the site of rapid loop production, so the 
many small cusps surrounding a crossing of the long distance parts of p± will 
actually end up on small loops, disconnected from the long string. Therefore, we 
expect a high tendency of the small loops to contain cusps. In loop trajectories 
with no kinks, the curves p± are smooth and must be periodic, of course. Hence, 
if they intersect this must happen an even number of times, and for the simplest 
motions this gives two cusps per loop oscillation [79] . However, we have seen that 
small loops should be dominated by large kinks and we expect them typically to 
possess one cusp per period. It might seem paradoxical that the cusp size (13.401) is 
parametrically larger than the loop size fl3.34l) . but this is another Lorentz effect: 
in the loop rest frame these sizes are of the same order. 

5.1.3 A model loop 

In this section we will display explicitly the trajectories of typical small loops 
as predicted by the arguments of Chapter 4, where loop formation was studied. 
Recall that the picture is that, due to the high fractal dimension of the curves p± 
on the unit sphere, large cusps break into many small cusps, each one having a 
tendency to form small loops around them. 

For such a "model" loop we take the functions p± to be short straight arcs on 
the unit sphere: 



These definitions hold for —f3 < u,v < (3, where 2f3 is the length of the loop. 
However, these functions can be extended to the whole {u, f )-plane by imposing 
2/3-periodicity. The jumps at u, v = P + 2Z/? correspond to kinks. The velocity of 
the cusp, V, and the period of the small loops, were determined in Chapter 3 
as average values over the ensemble of small loops. For a matter-dominated era 




(1, sm{Vu), 0, cos{Vu)) ~ (l, Vu, 0, 1 - V\y2) , 
(1, 0, sm{Vv), cosiVv)) ~ (l, 0, Vv, 1 - . 



(5.9) 
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these quantities are given by [59] 

V = 0.08{Gfi)-^-^H-^ , (5.10) 
f3 = 10{GfiY-h , (5.11) 

where t is the FRW time. 

The embedding of the loop worldsheet in target space is given by 

<A^^,v) = p^^{u)du + ^ j p^_{v)dv , (5.12) 

and is therefore aperiodic as a function of m or f separately: 

xl,Xu + 2(3,v) = x>:,,,Xu,v + 2f3) = xl,Xu,v) + {(3,0,0,sm{Vf3)/V) 

^ x>:,Au,v)+ (3 {1,0,0,1 -V^pye) . (5.13) 

However, from the above equation it follows immediately that x'^ ^ is indeed pe- 
riodic in a, with periodicity 2/3. 

One can use equations (15. 9p and fl5.12p to obtain snapshots of the loop as it 
evolves in time. For a fixed time t*, considered to belong to the interval [0,/3] 
for concreteness, the v coordinate may be written in terms of the u coordinate as 
Vu = 2t^ — u and then the string is parametrized by u alone. If m G [2t^, — P, P], 

Xw.s.(m) = ^ (~ cos(V^m), - cosiVvu), sin(V^M) + sin{Vvu)) ■ (5.14) 

On the other hand, if m G [P, 2t* + P] we define u'^ = u — 2P E [—P, 2t* — P] so 
that 

Xw.s.(m) = ^ (~ cos(\/0, - cos{Vvu), sm{Vu'J + sin(Vf„) + 2 sm{VP)) . 

(5.15) 

Figure [SH] shows several snapshots of the model loop (15. 9p . Strictly speaking, 
there is no self-intersection but at t = (or equivalently when u + v = 0) the 
string actually folds back on itself. The discontinuities of p± at odd multiples of 
P correspond to kinks on the loop traveling in opposite directions. When they 
meet the loop becomes degenerate and at the same time a cusp develops at the 
other end. However, adding a small (quadratic) perturbation generically renders 
it non-self-intersecting while the cusp is preserved. Thus, small cosmic loops are 
expected to be generally stable against fragmentation. 
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Figure 5.1: The left panel shows the xy-projection of four superposed snapshots of 
the model loop considered. For clarity a total translation by (1, 1) was performed 
and the x- and |/-axis are in units of [Vpy. The right panel gives a 3D view of 
the model loop at one of the two moments (the lightest curve in the left panel) at 
which the x?/-projection is degenerate. The 2;-axis was translated by —V(3 and is 
in units of (VPY. 



Small loops with sizes of order 2(3 have very large Lorentz factors. To show 
this, consider the impulsion of the loop, given by 



P^{u,v) = ^{pl{u)+if{v)) 



(5.16) 



and the center-of-mass velocity, which is obtained by averaging the spatial part 
of P^{u,v) over the worldsheet coordinates. 



dt f^^ da , , 



Therefore the boost factor is given by 



7 



Vi5 



(5.17) 



(5.18) 



For G[i = 10~^ the Lorentz factor fl5.18p is of order 10'^, as we have already seen. 
In its rest frame the nth harmonic will have frequency f^^^^ = j/S'^n, since (3 is 
its time periodicity in the FRW frame. Therefore, the observed frequency for a 
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loop moving at an angle 6 with respect to the line of sight is 



'rest 



n 



f, 



n 



7(l-|v|cos^) /3(l-(l 



n 



6 



) 



(5.19) 



cos6' 



In the special case of a loop whose motion is exactly aligned with the line of 
sight we see that the observed frequency is boosted by a factor proportional to 
7^. This effect can bring the lowest harmonics into the LIGO frequency band for 
Gn < 10~^, i.e. values of the string tension not yet ruled out by observations, and 
so the potential for detection arises. We shall consider such a possibility in the 
next chapter, but now let us address the scaling properties of the distribution of 
loops over their lengths. 

5.2 Scaling of the loop number density 

In this section we will study the evolution of the loop number density distri- 
bution. The equation governing the evolution of the number density of loops in 
a cosmic string network is a detailed balance determined by conservation of en- 
ergy. In what follows we solve this equation, first ignoring gravitational radiation 
and then taking into account this effect. But first let us make a few preliminary 
remarks and introduce some notation. 

5.2.1 Preliminaries 

Recall the many processes that can act to change the number density of loops 
with a given length /. Expansion of the universe causes the strings to stretch, but 
only on scales larger than the horizon size [65] . Therefore loops essentially do not 
grow in timeQ Self-intersections of cosmic strings can produce loops from long 
strings and also fragment loops into smaller ones. Of course, string intercommu- 
tation can lead to the absorption of loops back into the long string population 
but this process is strongly suppressed for small loops and so can be neglected. 
Finally, the coupling of matter in the form of cosmic strings to gravity means that 
the network radiates away part of its energy. This has been amply discussed in 
the literature [79l EHl EOl EH] and as a consequence loops shrink as time progresses 
and eventually disappear from the network. 

Obviously, the rate of loop production plays a crucial role in the detailed bal- 
ance: it persistently feeds the loop population. These get diluted by the expansion 

^ Loops bigger than the horizon distance are regarded as long strings. 
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of spacetime. On top of these effects, gravitational radiation continuously shifts 
the distribution toward smaller lengths as time elapses. In Chapter 4 we found a 
power law for the loop production function. The potential divergence at small dis- 
tances is supposedly cut off at the gravitational smoothing length scale. However, 
we will ignore the existence of this minimum length for the loops in what follows 
and see how far we can get. In addition, it has been suggested in (a fraction 
of) the existing literature that gravitational radiation is unnecessary to achieve 
scaling EZl- The results of the following section confirm this hypothesis. 

Let us now introduce the necessary notation for the rest of this chapter. Denote 
by dn{l,t) the number density of loops with size comprised between I and / + dl 
at a given time t. Then the quantity ^ has units of (length)"^. Thus, if it ever 
reaches a scaling regime during the cosmological evolution, it must eventually 
approach the following form: 

^ = t-'fm ■ (5.20) 

Borrowing notation from [38], define the length of loops in units of the horizon 
size, a = l/d\^. If the scale factor takes the form a{t) oc , the horizon size can be 
expressed as d-^^ = t/ {1 — u) and so the signature of a scaling regime in the loop 
number density is a solution of the form 

dfn ^ S{a) 
da a df^ 

In [38] S{a) was dubbed the scaling function. 

5.2.2 Evolution neglecting gravitational radiation 

As discussed in Section 14.31 if gravitational radiation is neglected the number 
of loops within a comoving volume changes only due to loop production [5^ : 



dt\ dl J -fHH^ \t 
Defining F{l,t) = ^{l,t) and inserting the power-law expression for the scale 



factor we obtain 



dl 

l+2x 

72/4 yt 



tF + 3uF = ^(U , (5.23) 
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which has the general solution 



FM : „,..j '7V + ^. (5.24) 



72(3z/ - 1 - 2x)t^ \t J f^" 

The function Qo depends only on the variable /. Noting that the inequality 1+2% < 
Su is satisfied both in the radiation- and matter-dominated eras, we conclude that 
at late times (t oo) the loop number density approaches 

dn , c f l^^^ ^ 



Thus, expressing everything in terms of a and (ih we indeed find that the loop 
number density approaches a scaling regime, i.e. it takes the form fl5.2ip with 

^ ^ 7'(3z/-l-2x) ^ ' 

This confirms that the loop number density does approach a scaling regime with- 
out taking into account gravitational radiation. The exponent 2% — 2 = —p takes 
the values —1.5 in the matter era and —1.8 in the radiation era. As we have 
already mentioned in Chapter 4, this is in good agreement with the numerical 
results of |38] who quote pmat = ^■4:1^q^j and prad = l-SO^g^g. However, these 
simple power-laws become good fits to the data only above a physical length £c, 
which was identified in [38] with the initial correlation length of the network. Fur- 
thermore, if the solution (15.261) were valid over the full range of loop lengths the 
total energy density would diverge in the UV since 

l^dl = d^^ J Sia)da . (5.27) 

We will now show that including the process of shrinkage of the loops due to 
gravitational radiation changes the power-law below the gravitational radiation 
scale, thus yielding a convergent total energy density in loops. 

5.2.3 Evolution including gravitational radiation 

Inclusion of gravitational radiation into the evolution equation for the loop 
number density introduces an extra term on the left-hand side of equation (15.231) 
because loops shrink at a rate given by (15. 3p . Therefore, we now have 



tF + 3uF - TGjjtF' = — - - . (5.2^ 
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Defining for convenience b = TGn, the solution of the above differential equation 
may be written as 

, , cb(l + bt)^^-* f bt V""'^" f bt 



72 \l + btj \l + bt 

Gil + bt) , , 

where B represents the Euler incomplete beta function and Q can be any function 
of the combination / + bt. 

By employing the Taylor expansion 

B(e;3u-1- 2x, 2x - 2) = _ ^ _ + 0(e) , (5.30) 

we find that the first term in fl5.29p behaves, for large I and fixed t, as ~ 
whereas the second term goes like ~ Q{^)- Requiring that the energy density in 
loops converges in the IR limit / — > 00 imposes that the general function Q{x) 
decays faster than a;"^. Therefore, using the expansion 

i? (1 - e; 3. - 1 - 2x, 2x - 2) = - " ^^^^-xy"'' . + o(-) , (5-31) 
^ ^ r(2x-l)r(3-2x) \ej ' ^ ^ 

the second term in fl5.29p is dominated by the first term as t — > 00: 
dn^^ , cbil + btf^-^ f bt V^'^^'f bt 

-(^ t) (rvn) b( —-■ ?>u-i-2x,2x-2 



dr'' 72 \l + btj \l + bt' 

(5.32) 

Once again converting to the variables a and d^, we obtain a solution of the 
form fICTD with 

(5.33) 

Now we can use the series expansions (15.301) and (I5.3ip to recover the limits 
for 'small' and 'large' loops. The separation between these two regimes is set 
by the gravitational radiation scale FG/x, and we find for a ^ FG/i the same 
result (I5.26p . whereas for a <^ FG/i 

7r(l-2x)csc(27rx) c 

'"^"^ F(2x)F(3-2x) 7^FG/.(l-z.)2+2x« ^"'^^^ 
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Since < 2% < 1 holds in both cosmological eras, the integral (15.271) is manifestly 
convergent and so the total energy density in loops is finite, as desired. 

Thus, taking into account gravitational radiation has the effect of bending the 
curves for aKTGfi so that the exponent in the scaling function becomes 2% — 1. 
Indeed, Figure 3 of that Ref. [38] does appear to show a certain range of the 
parameter a in which the scaling function behaves as such a power-law. At first 
sight this might seem intriguing since those simulations did not include the grav- 
itational radiation process directly. However, small loops behave like matter and 
so the expansion of the universe effectively shrinks the loops. Because the simula- 
tions only keep loops with sizes greater than a fixed fraction of the horizon, they 
are eventually removed from the game, hence emulating gravitational radiation. 
There, the minimum counting size was amin = 10"^ so we should expect the sub- 
gravitational radiation regime to set in for comparable scales or smaller. A more 
accurate estimate, equating the lifetime of the loops determined by gravitational 
decay to the lifetime set by the minimum counting size, yields abend = jE^c^min for 
the scale at which the bending would occur. Nonetheless, it is curious that this 
sub-gravitational radiation regime shows up also at early times in the simulations, 
when scaling is yet to be achieved. 

In conclusion, we have shown that under the very mild assumption that the 
energy density in cosmic string loops converges in the infra-red (IR), then the 
loop distribution evolves toward a scaling configuration in which it is also free 
from divergences in the UV as long as the process of loop decay by GW emission 
is considered. Note that the inclusion of gravitational radiation in our equations 
leads to a finite energy density of loops, even though the loop production function 
used as an input diverges at small scales. This just means that the rate at which 
loops are removed from the network is sufficiently high to balance the diverging 
loop formation. 
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Chapter 6 

Some implications for observation 



This chapter is concerned with observational signatures of cosmic strings. The 
tangible effects of these networks often depend on their detailed properties and 
the methods developed in previous chapters allow us to address some of these 
issues. 

We start by considering the implications of our model for the lensing by cosmic 
strings in Section 16. 1[ In this context we study both the typical distortion of 
the double images produced and the alignment of lensed pairs. The latter is 
related to the non-Gaussianity of the string distribution. We devote Section 16.21 
to the possibility of detecting (quasi-)periodic GW at the Laser Interferometer 
Gravitational Wave Observatory (LIGO) from the highly boosted population of 
small loops. 

This chapter is mostly drawn from Refs. [581 EH] , with Joseph Polchinski and 
Florian Dubath, respectively. 

6.1 Lensing 

Let us now consider the effect of the small scale structure on the images pro- 
duced by a cosmic string lens. A well known fact is that a straight string of 
tension /z has the sole effect of introducing a deficit angle in the transverse di- 
mensions, equal to SirGfj, [25], and this can induce identical images of sources in 
the background. However, once we include the possibility of small fluctuations 
propagating along the cosmic string the two copies will become distorted. 

Previous work has discussed possible dramatic effects [551 ES], including mul- 
tiple images and large distortions. We can anticipate that the rather smooth 
structure that we have found, which again we note is subject to our assumptions, 
will produce images with only small distortion. We will use our stretching model 
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Figure 6.1: Schematic representation of the system considered, with the string 
lens displayed along the Xi-axis and the distant source and the observer located 
at points S and O respectively. 



of the two-point function. If this proves incorrect one could apply the analysis 
using phenomenological values of x and A; for example, the extrapolation of the 
results of Ref. [37] give a smoother string, and even less distorted images. 

It is worth mentioning a recent occurrence that led to some excitement in the 
field. In 2003 a candidate for a cosmic string lens, CSL-1, was identified in the 
sky [97]. The two images of the putative lensed pair were identified with galaxies 
having apparently equal redshifts, morphologies and spectra. Furthermore, their 
angular separation placed the cosmic string tension at the higher end of the per- 
mitted range and several other potential candidates were identified nearby. Alas, 
a closer investigation with the Hubble Space Telescope revealed that the two im- 
ages indeed corresponded to distinct galaxies [98], showing no sign of a sharp 
discontinuity and the two images appearing distorted relative to each other. 



6.1.1 Distortion of images 



We quote the result of ref. [99] for the angular deflection of a light ray by a 
string, 



7±(y±) = 4:Gfi / da 



Vi/7(o)7(o) 



y± 



1 



Xll 



(6.1) 



t=toi<7) 



Here is the four-velocity of the unperturbed light ray, which we take to be 
(1, 0, 0, 1) as shown in Fig. 16. H is a reference point on this ray, and subscripts 
_L and || are with respect to the spatial direction of the ray. Also, x^{a,t) is the 
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string coordinate in terms of which 

F^-" = x^^x" - x^'x"" - ^r7^'^(x% - xf'x'p) , (6.2) 

and t^io) is defined by to(cr) = ^zi.'^-,^^ — Vz- 

In keeping with Section r2.3. II we separate the string locally into a straight part 
and a fiuctuation; we will keep the defiection only to first order in the fiuctua- 
tion. We consider here only the simplest geometry, in which the straight string is 
perpendicular to the light ray and at rest, so that P+ = — P_ = (1, 0, 0). To first 
order in the fiuctuation, x^{a,t) = (t,a,X2{<J,t),X3{a,t)). One then finds 



[a - yiy + 

72(y±) = -sgn{y2)A + AG fi da — — ^ , (6-3) 

J {(^-yi) +y2 

with A = AnGfi being half the deficit angle of the string. To the order that we 
work to is a constant, corresponding to the time when the light ray is perpendicular 
to the straight string. 

We can use our results for the small scale structure to calculate the two-point 
functions of the defiection. Let us focus on the local magnifications parallel and 
perpendicular to the string, given by basic lensing theory as 

AUV-A)97i. . 
M,(y,) = 1 — (y,), 

M.(y.) = (6.4) 

Do dy2 

Here Do and Di are the distances of the source from the observer and the lens re- 
spectively (these would be the angular diameter distances on cosmological scales). 
It is particularly interesting to consider the differential magnifications for the two 
images produced by a string, 

5Mi = Mi(y^) - Mi(yl) , dM^ = M^iy^) - M^iy'^) . (6.5) 

We take for simplicity y^ = (0, b) and y^ = (0, -b) for b = ADi{Do - Di)/Do\ 
this corresponds to the symmetric images of an object directly behind the string. 
Then 

5Mi = -5M2 = -— [ dax',{a,to)d, (^7^) ■ (6-6) 



^Thc region where the the Hght ray passes the string is small on a cosmological scale, so to 
use our earlier results we can locally set a = e = 1, dl = da, dt = dr- 
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From Section 12.3.11 we obtain 

(w+((T,t) ■w+(a',t)) = (w_(a,t)-w_(a',t)) =4(x^((T,t)x^((T',t)) 

= A[ia/tf^ + {a'/tf^ - {\a - a']/tf^^ . (6.7) 

Recall that Eq. (12.331) does not fully determine the two-point functions (16 ■7p . and 
in fact that the latter cannot be translation invariant. Here, we have fixed the 
ambiguity by defining the expectation value to vanish when a or a' vanishes (that 
is, at the point on the string nearest to the light ray); this amounts to a choice of 
how one sphts p-t into P-t and w±. 

As is stands, equation (16. 7p ignores the transition from a radiation dominated 
evolution to a matter dominated cosmology. In effect, if we take a small segment 
of string (in the matter era) and evolve it backward in time it will grow relative to 
the horizon size and eventually reach it, as we have alluded to before. Therefore, if 
the segment under consideration is short enough, it would have been comparable 
to the horizon distance before the radiation-to-matter transition, at a time teq, 
where the string dynamics is determined by the radiation era values. Hence, 
there is some critical length Ic and for \a — a'\ shorter than this one must take 
into account the transition. The details of this effect are collected in Appendix B. 
According to equation (lB.2p . the right-hand side of (16.71) just gets multiphed by 
a power of t/teq- From {x'^x'^ we obtaiij^ 

{ml) = {5Ml) = Xr(l - ^^'^) ^^(26/to)2x.(to/t^q)-2U-2x..U+2x. _ g) 
i-i COS! TT^'j- ) 

Plugging in the numeric values for x, A and ( and using a representative value 
for the dimensionless parameter Gfi ~ 10^'', we obtain a RMS differential magni- 
fication slightly below 1%: [fl 

{SM^y^^ ~ 0.009 . (6.9) 

Considering a smaller value of Gfi reduces the final answer. In any case, the take- 
away lesson is that the smoothness of the strings predicted by our model implies 
very small differential magnifications between the two images produced by a long 
cosmic string. 

^The lensing scale, which is provided by b, is typically of the order of 10~'^djy, thus much 
smaller than the critical length Ic- We can then safely use the formula for {x'2x'2) valid for the 
smallest scales over the full range of integration. Corrections from the longest scales will increase 
the final result, but not significantly. 

•^For this sample calculation the source was taken to be at a redshift of z ~ 2.5 and the 
string was placed midway between the observer and the source. The values to — lOGyr and 
icq — O.lOMyr were used for the intervening cosmological times. 
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6.1.2 Alignment of lenses and non-Gaussianity 

Another question related to short-distance structure is the ahgnment of lenses. 
Suppose we see a lens due to a long string, with a certain ahgnment. Where 
should we look for additional lens candidates? Previous discussions |100[ 1101] 
have considered the two extreme cases of a string that is nearly straight, and a 
string that is a random walk on short scales; clearly the networks that we are 
considering are very close to the first case. 

We keep the frame of the previous section, with the lens at the origin in and 
aligned along the xi-axis. Then as we move along the string, the RMS transverse 
deviation is 



The extension in the Xi direction is just /, so the RMS angular deviation is 



If we put in representative numbers, looking at an apparent separation on the 
scale of arc- minutes for a lens at a redshift of order 0.1, we obtain with the matter 
era parameters a deviation 6(f ~ 0.05 radians. That is, any additional lenses 
should be rather well aligned with the axis of the first. If the string is tilted by an 
angle ip to the line of sight, then projection effects increase 66 and 6(p by a factor 
1/cos-?/'. Of course, for lensing by a loop, the bending will be large at lengths 
comparable to the size of the loop. 

Lens alignment provides an interesting setting for discussing the non-Gaussia- 
nity of the structure on the string. If the fluctuations of were Gaussian, then the 
probability of finding a second lens at an angle 6ip to the axis of the first would be 
proportional to e"'''^^/^'''^ , and therefore very small at large angles. However, we 
have considered thus far a typical string segment, which undergoes only stretching. 
There will be a small fraction of segments that contain a large kink, and one might 
expect that it is these that dominate the tail of the distribution of bending angles. 

Let us work this out explicitly. Consider a left-moving segment of coordinate 
length a, and let Pla, r, k) dk be the probability that it contain a kink for which 
the discontinuity |p+ — p'_,_| lies between k and k + dk < k <2). There are two 
main contributions to the evolution of P. Intercommutations introduce kinks at a 
rate that we assume to scale, so that it is proportional to the world-sheet volume 




(6.10) 




~ v/^/4(x+l)(//t)^ = Sip . 



(6.11) 
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in horizon units, a? ear '^'^'^'^dr, and to some unknown function Q{k). Also, the 
expansion of the universe straightens the kink, k oc a~" [9]. TheiO 

^ = T^^'^'-''^"''^ag{k) + a^^{kP) . (6.12) 
OT a ok 

We set P to zero at the matching time tq defined by 

a = , (6.13) 

as in Eq. fl2.25p : earher kinks are treated as part of the typical distribution, while 
P identifies kinks that form later. For simplicity we assume that xq is small 
enough that the probability of more than one kink can be neglected. 
To solve this, define 

K = kT^', Q{a,T,K) = kP{a,T,K) . (6.14) 

Then 

TdrQ = ag{KT-'^')KT-^-''' . (6.15) 
This can now be integrated to give 

= ['^^k'"''3{k'). (6.16) 



Here ko = k{xo/a)^'^^^^^''"°^'^ and a = cx/t^^^'''^\ Note that a is just a constant 
times l/t, so the probability distribution scales. The source g{k) vanishes by 
definition for A; > 2, so /cq > 2 is equivalent to fco = 2. 

Rather than the angle 6(p between the axis of the first lens and the position 
of the second lens, it is slightly simpler to consider the angle 56 between the two 
axes. In the small fluctuation approximation this is just X2{(t), and so the RMS 
fluctuation is 

59= ^/A/2{l/t)'' . (6.17) 

In the Gaussian approximation, the probability distribution is ^^^^^ 1'^^^ and so is 
very small for large angles. On the other hand, a large angle might also arise from 
a segment that happens to contain a single recent kink. Treating the segment 



^ Other effects are often included in the discussion of kink density, such as the removal of 
kinky regions by loop formation [57 , but these have a small effect. 
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as straight on each side of the kink, the probabihty density is then precisely the 
function P{a, r, k) dk just obtained, with k = 66. If we consider angles that are 
large compared to 69 but still small compared to 1, the range of integration in 
the solution fl6.16p extends essentially to the full range to 2 and so the integral 
gives a constant. Then 

P((T, r, k) dk oc k-^/^dk = 6e-^'^d6e . (6.18) 

Thus the tail of the distribution is not Gaussian but a power law, dominated by 
segments with a 'recent' kink. One finds the same for the distribution of 6>^. It is 
notable however that the exponent in the distribution is rather large, roughly 4 in 
the matter era and 10 in the radiation era. Thus the earlier conclusion that the 
string is rather straight still holds. The sharp falloff of the distribution suggests 
that a Gaussian model might be a good approximation. 

6.2 Periodic gravitational waves from small cos- 
mic string loops 

We have seen in the previous chapter that the small loops which our stretching 
model predicts to be abundant are characterized by very high velocities. Such a 
strong Lorentz boost can push the natural oscillation frequency of the loops into 
the frequency band of GW interferometers. The small loops may then be potential 
sources for these experiments and if so they should give rise to somewhat special 
GW signatures. In this section we study this possibility. 

One of the most promising ways to test the presence of a cosmic string network 
in our Universe is through its gravitational wave (GW) emission [88l |92| [93| [102j . 
One can expect a GW background from the network [103[ I104[ [9l] and on top of 
this the loops typically form cusps which emit strong bursts of radiation [271 180] . 
The presence of a population of highly boosted loops leads to a peculiar GW 
signature and therefore deserves a specific study as these may be potential sources 
for present and planned gravitational wave detectors. The large Lorentz factor 
has some consequences on the analysis of GW detection which were not previously 
considered in [Ml EH 1021 El ED] . 

All cosmic string loops are expected to have (quasi-) periodic behavior, but 
one of the key features of small high-velocity loops is the fact that their observed 
period may be as short as a few mili-seconds. One may even expect that the first 
few harmonics of the loop enter the frequency band of the GW interferometers. 
Such a case is very promising since, on one hand, the lowest harmonics emit 
the strongest signal and furthermore it has little dependence on the exact loop 



77 



CHAPTER 6. SOME IMPLICATIONS FOR OBSERVATION 



trajectory: we can expect more robust wave-forms than those originated from 
the single cusp case (whose behavior under back-reaction or in the presence of 
small-scale structure has been questioned |105l 1671 1761 1821 [106]). On the other 
hand, if this frequency is short enough in order to enter the GW interferometer 
band, this allows to search for periodic GW signals accumulating a large number 
of loop periods. As a result the strain sensitivity will be improved by a huge factor 
corresponding to the square root of the ratio between the observation time and 
the lifetime of a single cusp. 

The fraction of the small loop population moving directly toward us sees its 
GW strength boosted. However, the large Lorentz factor also pushes most of the 
emission of gravitational radiation into the direction of the loop motion. The fact 
that the signal is potentially strong enough is not sufficient to guaranty observa- 
tions; we need the fraction of sources meeting the observability requirements to 
be numerous enough in order to get a reasonable probability for detection. 



We now compute the GW spectrum emitted by a typical small loop. The 
model loop described in Section [5. 1.31 will be considered as a representative of this 
class. To distinguish between the location of the cosmic string and a generic point 
in spacetime we shall use coordinates ^ for the former and coordinates for 
the latter. 

A cosmic string acts as a source term for the gravitational field through its 
energy-momentum tensor. For a classical string it is given by [6] 



with X = (t,x), while the four- vector x„.s.('U,f) gives the spacetime location of 
the world-sheet point with coordinates {u,v). Since the loop has 2/3-periodicity 
we choose the world-sheet to be a strip {u,v) G {—oo, oo) x [—/?,/?]. 

Furthermore, for a source of size ~ d localized around the origin, the trace- 
reversed metric perturbation in the local wave zone (r = |x| >> d ) is given, in 



6.2.1 Emission of GW 



T^'^(x, t)= fi {p1p''_^ + p^^pI) 5^ {x - v)) dudv 



(6.19) 



the time domain, by P 11071 1108] 





[p+p"- + P-P\) 6 (t - r + Xw.s.(m, f ) ■ n - a;° .sX"", v)) dudv , (6.20) 
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where n = x/r. When analyzing high frequencies originated by cusps and kinks 
one must be careful with the fact that the dominant contributions to p'^ in the se- 
ries expansion is a pure gauge term, as was first noted in [80] . However, converting 
to transverse traceless (TT) gauge eliminates this termo In any case, we will be 
interested in the lowest harmonics since these are the frequencies which contribute 
the most to the observation rate. Performing a temporal Fourier transform one 
obtains 

%'"'{iu, x) = ^e'^" j {j/Ip''_ + e'"^'^-""- rfwrft; . (6.21) 

Now, decomposing the integral over u into 

/oo p{2m+l)l3 
du = y2 ^« ' (6-22) 

■oo ^.^-7, J(2m-l)0 



using property (15.131) and defining Q as the angle subtended between n and the 
2;-axis (along which the loop is traveling) we get 

J -13 J -13 

V ' 

jMf (6I,</,) 

= |v| cosg) 2^ n /3 (1 - |v| cosg) ) ^ 

= ^e^"7iWE'^(^-2^/"(^))^^'^(^''/') ■ (6-23) 

n 

As a check, note that we have recovered the discrete set of frequencies obtained 
in f l5.19p . The computation of I^'^{u,v) involves the following integrals: 



u,3 



jn^ ^ l"^ ^j^2niMe)(^^{l-cose)-smecoscP^+cose^)^^ ^ ^g 24) 

-13 

4", = ^ie2.^/"W(f (l-cos9)-sin9.in0^+cos(?^)^^ , J = 0, 1, 2 . (6.25) 

/3 



^One can also eliminate the gauge term explicitly by replacing p^. by Pj. = — (1, — n). 
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Before we move on to computing fl6.23p one more consideration is in order. The 
interaction between a GW and a detector is usually described in the TT gauge 
and it is convenient to rotate the coordinate frame in order to match the observer 
description in which the GW arrives along the 2;-axis, so that i"'^™ = n. The 
rotation matrix that converts between the source frame and the observer frame is 
given by 

(cos 9 cos (j) cos 9 sin cj) — sin 9 \ 
-sin0 COS0 . (6.26) 

sin 9 cos (p sin 9 sin (p cos 9 J 

Rewriting h^^, in the TT gauge is performed by means of the projector Ajj fc/ 
defined by (see section 10.4.15 in |108] ) 

^ij,kl = PikPjl — -PijPkl , (6.27) 

where P^ = 6ij — riiUj. Then, 

hj^ = KjMRkk'Rwh'v = \ I , (6.28) 

and we can define the observed strength as 



h= ^\h+\^ + \h^\^ . (6.29) 

Note that the main dependency of h on the direction of observation comes from 
the fi{9) factor in Eq. f l6.23p which, for small values of 9, causes an enhancement 
by a factor (V 13)^"^ (~ 10^ for G^i = 10^^ ), as was discussed in Section [5. 1.31 The 
remaining angular dependence arising from I^'^{9,(l)), for the first harmonic, is 
shown in Figures [6l2] and [6731 One can observe that this gives some enhancement]^ 
for 9 ~ 7r/2. However, this corresponds to the range in the parameter 9 where 
the large Lorentz boost enhancement is lost. Let us anticipate that in order to be 
able to detect the loop we need a big enhancement of the emitted GW (in other 
terms, we are able to 'see' only loops moving nearly in our direction). Therefore 
any sub-leading effect at large 9 may be neglected and in the remainder of this 



^This enhancement is a consequence of the rotation between the source frame and the observer 
frame. This is more clearly seen by taking = 0, in which case the dominant contribution to 
h/ fi{6) is proportional to sm^{6). 
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Opi g J, Opi 



Figure 6.2: The angular dependence of the first harmonic (n = 1). The plots 
show the quantity h ■ 2^;^^^—^t^^- We have set Gjj, = 10^^. The right panel is a 
closeup on the value of 9 of order V (3 ~ 1.5 ■ 10~^. 



study we will discard this sub-leading dependence. This leads to the following 
form (omitting the higher harmonics) for h: 



1.3- 



r (1 

TT 



V COS! 



xe 



dx 



■0.5 



r 



cos 9 + 



^cos9 



(6.30) 



where the final factor of 0.5 shifts the value of the integral at ^ = to its actual 
minimal value (see Figure in order to obtain a lower estimate for the strain h. 



6.2.2 Visibility of the small loops 

Having computed the GW spectrum, we now want to determine whether or 
not current or planed GW detectors can observe the small relativistic loops under 
consideration and, if so, for what range of the parameters. This will depend 
on the cosmology and therefore we must fix the dependence of t and r on the 
redshift z. To this end, we consider a flat universe with the cosmological constant 
contribution included in the cosmology. 

Any GW detector is characterized by its frequency window [/_;/+] as well 
as the minimal amplitude needed for a detection /idet(/)- For LIGO's second 
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Figure 6.3: Same as Figure ESI but after averaging over 0. Note that for VP < 
6 < 71 — V(3 the curve is approximately proportional to sin^ 9. This dependence 
arises from the rotation and posterior conversion to TT gauge. 



science run (LIGO S2) the frequency window is [160; 728.8]if2;. Furthermore, the 
sensitivity to a continuous signal from a 10-hour search can be modeled by 



where q ~ 0.6 and /iieo ~ 10~^^, see |109j . LIGO is expected to be 10 times more 
sensitive and Advanced LIGO to bring another factor of 10 (see for example [llOj ). 

In order to detect a given loop there are 3 conditions which must be met: 
first, the frequency of the incoming GW produced by the loop has to lie within 
the detector window; second, its amplitude has to be above the limit /idet(/); 
finally, the observed lifetime of the loop, Tobs, must be large compared to the 
time span of the experiment, Tobs, to insure that its periodicity does not change 
significantly during observation, in which case it would quickly drop out of the 
frequency window anyway. 

Inspection of equation fl5.19p shows that the frequency of the received GW 
depends on the angle 6 subtended between the line of sight and the velocity of 
the loop. Thus, the strategy is to first fix Gjj, and determine the range of angles 
under which observation is possible, for a given z. We shall give the different 
bounds in terms of the cosine of the corresponding angle. Then, for each value of 
G/i, a combination of the above constraints determines an observational window. 
An example with G/i = 10~^ is shown in Figure 16.41 Let us point out that, if 




(6.31) 
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z 



Figure 6.4: Observational window for LIGO and for G/i = 10^^. We plot cos(6') 
as a function of the red-shift. The solid curves represent the frequency window, 
the dashed curve is the sensitivity bound and the dash-dotted curve corresponds 
to the lifetime constraint. Note that, for such a value of Gfi, we can observe loops 
up to 2; ~ 9 when the direction of the motion of the loop lies on a cone of aperture 
~ 2 ■ 10~^rad ~ 0.7' around the line of sight. 



our assumptions are correct, one can potentially observe these loops at LIGO if 
10"^° < Gfi < 7 X 10"^. Finally, the integration over the redshift can be performed 
to determine the expected number of detections. The capacity of observing such 
loops (or placing constraints on the string tension) will then depend on their num- 
ber density in the universe, and this quantity will be worked out in Section 16.2.31 
But before we do so, let us turn to the analysis of the several constraints mentioned 
above. 



Constraint from the frequency window 

In order to determine the observed frequency of a GW produced by a small 
loop we must correct equation (15.191) for the cosmological redshift, thus obtaining 

/niobserved = ^(1 + ^) (l - COS 6 {1 - (3^ 6)) ' ^^'^^^ 
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The bounds /± on the observed frequency lead to corresponding constraints on 
the angle 6 which can be expressed as cos^^_ < cos^^ < cos^^+, with 

1 - 



cos^± = (6.33) 

Constraint from the sensitivity 

The amplitude h of the GW wave produced by a small cosmic string loop was 
obtained in (16.301) . Taking into account the cosmology and converting to the time 
domain, the strain caused by the lowest harmonic takes the following form: 

Ht, e) = !,f,^C^L,, , , cos(27rA(e)t - ^o) , (6.34) 

(1 — cos 6^(1 — V^p''/6)) r[z) 

where ipo is an (arbitrary) phase. The condition to be satisfied in order to be able 
to detect the GW is then h(t, 6) > /idet(/i| observed (6')). Solving for the angle yields 
a lower bound, cos 6* > cosOh, where 

1.3 GfiV'^l3^+''{l+z)i(imHz)" ^ 



COS 9h = ^ _ iy2p2 • (6-35) 

Constraint from the lifetime 

Another feature that must be taken into account is that cosmic string loops 
shrink with time due to the loss of energy in the form of gravitational waves 
themselves. We have already seen in Section 15.1.11 that the lifetime of a loop of 
size 2(3 in the FRW frame is r = (2/5)(rG'/i)~^. Since the loops are moving with 
a relative velocity which makes an angle 6 with the line of sight their observed 
lifetime, including the cosmological redshift, is 

Tabs = (1 - |v| cos^^)(l + z)t = (1 - |v| COS0)(1 + . (6.36) 

Given the extremely high velocities of the loops, one finds that for very small 
angles 6 this apparent lifetime can be very short since it is suppressed by a factor 
of 7^. In that case, if we want to consider loops whose frequency does not change 
significantly during an observation time Tobs we need to guaranty that robs > Tohs, 
which translates into the constraint cos 6 < cos 6^ , where 
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This bound has the same form as the bound fl6.33p coming from the frequency 
window. The hfetime constraint can thus be recast as a constraint on the frequency 
window instead. Combining equation (16.331) and (16.371) we obtain an effective 
maximal frequency given by = minj/^, //,}, with 

2n 

h = ^ . (6.38) 

i (jr /ii obs 

The hfetime constraint has the effect of introducing a G/x-dependent maximal 
frequency. For the first harmonic {n = 1) and a time span of Tobs = 10 hours in 
the settings of LIGO, this does not affect the original frequency window for Gfi < 
1.5 X 10~^. Beyond that, the maximal frequency is given by fi, which decreases as 
Gfi is increased, until it finally closes the frequency window for Gfi ~ 6.9 x 10~^. 
If we want to observe or constrain higher values of G/i we have to look at higher 
harmonics or make a shorter run. Both possibilities help relax the bound coming 
from the lifetime but on the other hand they also tighten the sensitivity constraint. 

6.2.3 Loop number density and expected number of de- 
tections 

For the purposes of computing the expected detection rate we shall assume 
that the string network is in a scaling regime. The rate at which the long string 
is converted into small loops can be obtained from equation (I4.27p . Inserting the 
numbers for a matter-dominated era, where the energy density in long strings is 
Poo ~ 4/xt~^, one finds that 

di^\ 4 Vol 

(6.39) 



where Vol is the Hubble volume. We obtain the number density in small loops 
by multiplying the above equation by the lifetime and dividing by their length 
as well as the volume factor. Since the lifetime of the loop fl6.36p introduces a 
dependency on the angle of observation 6, the apparent number density is given 
by 

n{9) = (1 - |v| cos^)n^w = (1 - |v| cos^^)^-^^ . (6.40) 

5 I Gp t'^ 

Now we can include the constraints from the Section 16.2.21 Only loops moving 
along certain directions and within certain distances emit GW which enter the 
frequency window with an amplitude high enough to be detected. The range of 6 
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for which the detector can 'see' the loop is given by cos 6'mm < cos6' < cos 6'maa;, 
where 

cos 9max = niax { — 1, min {cos ^^_|_, COS 6*^}} , 

cos = min {1, max {— 1, cos6'_, cos6'/j}} . (6.41) 

Therefore, averaging over the sphere gives 

- / n{e)smede = 

2 If) 

= ^ (^COS 9max " COS 6'„i„ - -y- COS^ 9^ax + ^ COS^ 9min^ ripRW ■ (6-42) 

V ' 

Finally, to compute the expected number of observations we must account for 
the non-trivial evolution of the Universe. To this end we will follow Ref. [93j, 
where the cosmology is encoded in the function h{z) = H{z)/Hq which expresses 
the time-dependent Hubble factor in terms of the redshift: 

h{z) = ^VL^{1 + zf + VLril + zY + . (6.43) 

The coefficients fi^, and are the present energy density contributions from 
matter, radiation and cosmological constant, respectively. These add up to unity 
in a flat universe and Vim = 0.25, Vt^ = 4.6 x 10~^. The comoving variables are 
expressed as t = HQ^ipt{z) and r = HQ-^ipr^z), where 

..W = f^. (6.45) 
Then, the number of observed loops is simply 

N = I C{0^ax,9min)nFRwdV , (6.46) 



where the comoving volume is given by 



dV=^—^^dz. (6.47) 
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Figure 6.5: Expected number of observations for a run of Tobs = 10 hours. The 
dotted, dashed and continuous curves represent LIGO S2, LIGO and Advanced 
LIGO, respectively. 



and Ho represents the current value of the Hubble parameter. Therefore, we 
obtain the following expression for the expected number of observations: 



As a result, we report in Figure [631 the expected number of detections for a run 
of 10 hours with the different versions of LIGO. We can understand heuristically 
some qualitative features of the curves. The rise of the curves for increasing Gfi 
comes mainly from the dependence of the GW strength on this parameter (see 
equation (16.301) ) which allows to observe a larger volume0 while the abrupt cutoff 
at large values of Gfi has its origin in the closing of the frequency window by 
the upper bound from the lifetime constraint. Another effect that contributes to 
this fall-off is the fact that the regions allowed by the frequency window and the 
sensitivity bound do not overlap anymore for high Gfi. 

Note that reducing the observation time has 2 effects: the first one is to drop 
the GW strength relative to the detector sensitivity since /iigo is proportional to 
T~^^J'^ |109j . On the other hand, this allows the detection of short-lived loops at 

''Equation (|6.35p shows how the sensitivity constraint becomes less restraining with increasing 




(6.48) 



87 



CHAPTER 6. SOME IMPLICATIONS FOR OBSERVATION 



higher redshifts by relaxing the hfetime constraint and consequently we can scan 
larger values of Gfi. However, this leaves the expected number of observations far 
from being of order 1 or larger. 

Note also that we have only taken into account the first harmonic. Adding 
higher harmonics has the effect of expanding the curves in Figure 16.51 to higher 
values of Gfi but this does not cause the curves to rise. For specific combinations 
of the parameters Gfi and z it is possible to observe more than one harmonic 
in the detector frequency window. This may provide (through a specific search 
for direction correlated signals) better sensitivity. However, this is certainly not 
enough to get an expected number of observations of order 1. 

6.2.4 Summary and discussion 

In conclusion, we have shown that the observation of boosted cosmic string 
loops is, in principle, possible for very interesting values of Gfi in the range 
IQ-^o — 10^®. However, the chances of such an observation during a 10 hour 
period are slim due to the low effective number density of such loops. Even 
though a large boost factor can bring the typical frequency of the GW produced 
into the LIGO frequency band, this will only occur if the loop in question is mov- 
ing along a trajectory very close to the line of sight {6 < {Vf3)'^), which effectively 
cuts down the apparent number density of loops with these characteristics. An 
improvement can be obtained by increasing the duration of observation, thereby 
gaining in sensitivity of the detector. If one keeps increasing the observation time, 
the apparent lifetime of the loops will eventually be reached. We point out that, 
due to the production mechanism of the small loops (most of the loops that are 
produced around a large cusp move in coincident directions), we do not expect a 
Poisson distribution. As a result, while longer times of observation lead to higher 
probabilities of detection, the expected number of loops observed should grow 
less than linearly with Tobs- 

We should stress that the possibility of observing GW from cosmic strings we 
have considered concerns only the continuous emission of GW by the lower har- 
monics of the loops. This is in contrast with [92l [93l [271 EO] which focus on the 
burst of GW produced by the cusps and kinks. The results we have presented here 
depend both on the exact characteristics of the population of small loops (namely, 
their size) and also on the type of cosmic string network. We have assumed the 
distribution of loop size to be sharply peaked at the gravitational radiation scale 
but a decaying power law would reproduce a more realistic population. Also, 
cosmic superstrings can have reconnection probabilities as small as P ~ 10~'^ [71] 
and the number density of loops is proportional to a negative power of P. Nu- 
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merical studies I112j suggest values between —1 and —0.6 for the exponent. 
This would lead to an enhancement of the probability of detection but not quite 
enough to obtain iV ~ 1. On the other hand, our conclusions are independent of 
the presence or absence of high frequency features on the loops, i.e. cusps and 
kinks. 

Note also that the discrimination between a signal from a boosted loop and 
other periodic signals is not difficult since the former is expected to have a rising 
frequency and a decreasing amplitude. This is in contrast with both spinning stars, 
which are expected to spin down, and with mergers, which see their GW amplitude 
growing with time. However, a complete understanding of the waveform is not 
possible without a description of the higher harmonics. Obtaining a template 
for such waveforms would be a valuable achievement which would allow direct 
searches of GW emission from cosmic string loops. 

Finally, we have restricted our analysis to the LIGO frequency band. For the 
lower frequencies of the LISA band the potential sources are loops with high Gn 
at very large z. However, in order to get a sufficient signal-to-noise ratio we have 
to consider longer observation times and then the corresponding lifetime bound 
closes the observational window. 
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The evolution of the interest in the subject of cosmic strings has shifted since 
their original proposal. In the early eighties cosmic strings were regarded as a 
possible explanation for structure formation in the Universe, but more recently the 
interest has been revived due to the connection with String Theory. In particular, 
such objects could provide an observational window into this promising unifying 
framework. If this is the case and cosmic strings are ever detected, detailed 
knowledge of their networks is crucial to extract the real physics contained in the 
data. The work presented in this dissertation represents a step toward determining 
properties and effects of cosmic strings. 

Over 25 years of investigations have not led yet to a complete understanding 
of these networks. This is a clear indication of their high degree of complexity: 
the non-linearity of the equations of motion, the presence of a large hierarchy of 
scales and the interplay between the several mechanisms that take place. The 
issue of the build up of small scale structure on such systems stands out among 
these uncertainties. 

Thus, our approach has been to focus on a range of scales for which one mecha- 
nism (stretching due to the expansion of the Universe) becomes dominant and the 
equations of motion become approximately linear. The analytic model presented 
here therefore addresses only the short distance properties of the networks, avoid- 
ing the complex dynamics at scales comparable to the horizon size. The numerical 
simulations are needed to address this highly non-linear evolution at large scales 
but analytic methods are important to separate real effects from transients and 
numerical inaccuracies. Purely analytic or numerical treatments are frustrated by 
the complexity of these systems but the conjunction of both methods has led to 
some advances, namely on our understanding of the distribution of cosmic string 
loops. 
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How successful have we been with our efforts? Using the stretching model 
we have been able to fix the functional form of the two-point function between 
points on a long string, confirming their smoothness at short distances. The ex- 
pression we have found fits the data from numerical simulations on a range of 
length-scales below the horizon but there is some discrepancy at smaller scales 
which still begs for a unanimous explanation. The model developed also enables 
a more accurate determination of the scale of string smoothing by gravitational 
radiation. This is parametrically smaller than the horizon size by a power of the 
dimensionless string tension, yielding a large hierarchy of scales. The exponent 
is related to the two-point function. Furthermore, a novel picture for the distri- 
bution of loop sizes has arisen from both numerical and analytic investigations. 
Traditionally, loops were considered to be produced at a single length-scale. More 
realistically there is a power law distribution accumulating loops at the gravita- 
tional radiation smoothing scale and a second peak roughly an order of magnitude 
below the horizon distance. Unfortunately, our methods have not allowed us to 
determine the normalization of the distribution to confront with the simulations. 
Furthermore, the issue of loop fragmentation still remains a little obscure. We 
have also demonstrated that the loop distribution approaches a scaling regime in 
matter- and radiation-dominated eras, a point closely related to the long-standing 
question of the presence or not of scahng in the small scale structure. 

More generally, our work has unveiled several correlations among the two-point 
function, loop production, fragmentation, and loop velocity. These connections 
have led us to examine some observational consequences of cosmic string networks. 
Regarding the lensing phenomenon by cosmic strings we now have a better sense 
for the expected magnitude of the distortion of lensed pairs. On the other hand, 
detection at LIGO of gravitational waves generated from the low harmonics of 
small cosmic string loops is an exciting possibility for an interesting range of the 
string tension. 

In conclusion, despite the simplicity of the equations of motion, the evolution of 
cosmic string networks turns out to be inherently (and persistently!) complicated. 
Our efforts have been directed toward the study of the simplest possible kind of 
network and even so we still fall short of a complete understanding. Nonetheless, 
with the stretching model we have made some progress in relating properties of 
the system at different scales and in identifying subtle relations between disparate 
quantities characterizing the network. 

Gradually we are evolving toward a clearer understanding of the properties of 
cosmic string networks. Hopefully, the desired level of knowledge is an attractor. 
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In this appendix we include the calculations that support the estimation of 
the second term in the expectation value fl4.20p obtained in Section 14.2.21 namely 

((5(L;-Ll)|detJ|)pX^^ . (A.l) 

We shall assume that the functional probability distribution for the fluctuations 
is Gaussian with variance given by the two-point function. Furthermore, we see 
no reason for a strong correlation between the (5-function and Jacobian factors in 
the expression abovqj, so we take the product of their averages. 
For the 5-function, a Gaussian average gives 



oo 

1 



^ ^,({L^-Li))^-^^((L;-Ll)2)^-^^^((L;-Li)3)^+. 

27r 

27T 



(A.2) 



The subscripts c in the second line refer to the connected expectation values, 
obtained by contracting the Gaussian fields uj± with the propagator (14.111) . We 



^The whole point of the separation of variables in Section 14.2.21 was to accomplish this 
factorization. 
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have defined 

Rix) = 
c{x) = 



u+l/2 pu+l/2 /2+4x 

du' / du" (cu+K) • u:^{u")f = - ,(A.3) 

u-l/2 Ju-l/2 1^ 



_ l+2x I 1 
4(l+2x)2 V l+4x (1+X)2 



^ 4r^(2+2x) 

3+4x r(4+4x) 



(A.4) 



Numerically C{x) = (0.0049,0.0106,0.0111) for x = (0.1,0.25,0.5). 

We now consider the expectation value of the Jacobian, {| det J|)pX^^. The 
matrix itself can be written more concisely if we shift the coordinates u and v by 
1/2. In that case the quantities L± depend on / only through one of the integration 
limits. Consequently, the third row (obtained by differentiating L_|_ — L_ with 
respect to /) simplifies and we arrive at 



cj^(/) -cj^(O) 



(A.5) 



For notational simplicity we have set m = and f = / at the end. 

In the Gaussian approximation the relevant probability distribution is 



p±(a;±(/),a;±(0)) = — — ^exp 



--(B^)^M±B'i 



(A.6) 



where the index i is summed over the two coordinates x and y since the fluctuations 
obey the orthogonality condition z ■ uj± = 0, and 



=(0 

:(0) 



(A.7) 



(That is, the columns and rows of the 2x2 matrix M correspond to the points 
and /, not the index i.) As usual, the whole distribution is determined solely by 
the two-point functions which have already been determined in (14.111) and (I4.14p . 
Thus, one finds that 



2m^f 1 X 

A{l-x)\x 1 



detM 



± 



VJ i-x 

We now have all we need to write down an expression for ( | det J | ) . Before 
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we do so, let us perform a simplifying change of variables: 



(A.9) 



under which the expectation value of the loop Jacobian takes a compact form: 



F(X) + 



X 



1 + X 



G(X) 



27r4 Vt, 

Here we have defined an 8-dimensional vector 

X = (x+, x_, r+, z+, w_) 

and two functions of it: 

F(X) = (VVW^- - Y_W+) {Y^ + Wl-Y^~ Wl) 
G'(X) = 2{W+W^ -Y+Y^){X+- X_){Z+- Z_) - 
+ {Y+W- + Y^W+) [(X+ - X_)2 - (Z+ - 



(A.IO) 



(A.ll) 



(A.12) 



It is challenging to proceed further analytically but one can find upper and lower 
bounds for the value of the integral in (lA.lOp by noting that 



\A\ 



\B\ 



< 



\A + B\< 



\A\ 



\B\ 



(A.13) 



Indeed we can compute analytically the integrals of the functions F and G sepa- 
rately. Here we show how to accomplish this for the integral of -F(X). Under the 
following change of variables 



(A.14) 





= r cos 6 




= r sin 6 


Y_ 


= pcosip 




= psimp 
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that integral becomes 



/ 



27r 



TT / d9dip\ sin((y9 



drdp r p \p — r | e 



r —p 



TT 



de2 dtp sm{(p -9) - = 2tt^ . (A. 15) 

Jo Je 4 



A similar (nevertheless, more cumbersome) game can be played with the integral 
of G(X). It involves more complicated changes of variables but the result is just 
as simple: 

y"c/«Xe-^' |G| = TT^ (A.16) 
Thus, the expectation value of the Jacobian is bounded by 



1 i-x 



TT 2 + 2x 
Numerically, we find that 

^-'(t/0'^(|detJ|) 



<^-'(t//)^^(|detJ|) < 



1 , i-x 



TT 2 + 2x 



0.57 during the radiation era 
0.48 during the matter era 



Finally, we can combine results (1A.2I) and ( 1A.18I) : 



(5(L;-Ll)|detJ|)pj 



t \t 



(A.17) 



(A.18) 



(A.19) 



where rj ^ 2.3 in the radiation era and rj ^ 1.3 in the matter era. 
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Appendix B 

The matter-radiation transition 



In this subsection we will assume that our stretching model is actually valid 
down to the scale where gravitational radiation sets in. If loop production or other 
relatively rapid processes are actually determining the small scale structure then 
this appendix is moot. 

We have noted that at very short scales we see structure that actually emerged 
from the horizon dynamics during the radiation era. Thus we should take the 
radiation-to-matter transition into account in our calculation of the small-scale 
structure. At the time of equal matter and radiation densities, 

1 - (p+((T, Teq) ■ P+((t', T,^)) ^ A(^eqAeq)'^^ , (B.l) 

where leq is the length of the segment between a and a' at teq- Assuming that 
the transition from radiation-dominated to matter-dominated behavior is sharp 
(which is certainly an oversimplification), we evolve forward to today using the 
result (12.231) . The right-hand side of Eq. (]B.1|) is then multiplied by a factor 
(V^)"^''""'"- III terms of the length today, / = /eq(V'^eq)^", we have 

1 - (p+(a, r) ■ p+(a', r)) ^ Ml/tY>^^{t/t,,)-'^--'^^<-^'>^'- . (B.2) 

This expression applies to scales l(t) that, evolved backward in time, reached 
the horizon scale dn before the transition occurred, i.e. at a time t^, defined by 
l(t^) ~ dfj such that < teq- For longer scales, which formed during the matter 
era (for which > teq), we have simply 

1 - (p+(^, r) ■ p+(a', r)) ^ AUl/t)'''- ■ (B.3) 

The transition between the two forms occurs along the curve determined 
by the intersection of the two surfaces (IB.2P and flB.31) . This determines the 
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0.01 - 



^ 0.006 - 



V 0.004 - 



Ol ] ] ] ] ] ] ] ] ] 1 

0123456789 10 



Figure B.l: Structure on the string, as a function of the length / at present 

time (sohd curve). On scales larger than the critical length /c ~ 3 x lO'^dn 
the structure is determined by the matter era expression. On scales below Ic 
the transition result flB.2l) gives an enhanced effect. The dashed curves show the 
extrapolations of the two relevant expressions: on small scales the actual structure 
is enhanced relative to the pure matter era result. 



critical length at the time of equal matter and radiation densities, lc{teq) = 
{At/ AraY^^'^^"'~'^^''^teq- In terms of the length at some later time t, the transi- 
tion occurs at 



SO that the transition scale at the present time is /c ~ 3 x 10^^ dn (Fig. IB. 11) . The 
transition result (1B.2P implies more structure at the smallest scales than would 
be obtained from the matter era result, by a factor (IJl)'^^^-^-^'^^ ~ (/c/0°'^- 

Of course, precise studies of the small scale structure must include also the 
effect of the recent transition to vacuum domination; this period has been brief 
so the effect should be rather small. 
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